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Summary  " 

Card  programmable  calculators  together  with  suitable  prop 
grams  can  easily  provide  numerical  answers  which"  formerly 
took  large  computers  and  volumes  of  tabular  or  plotted  outputs 
for  everyday  reference . Presented  here  are  derivaUons  and  , 
program  listings  for  th-:  HP-65  and  HP-67  calculators  to  pro- 
vide the  User  with  commonly  used  radar  detection  performance 
data,  ' / ' - ' 

, *,  i , 

Both  fixed  threshold  and  adaptive  threshold.  C FAR  detection 
with  noncoherent  integration  are, covered.  Recursive  programs 
for, the  general  chi-squared  target  fluctuation  distribution  are 
treated  as  well  as  faster,  running  recursive  programs  for  the 
Swerling  target  Cases  I r.  IV.  Also  included  are'  fast  programs 
for,  the  required  detection  signal-to-ripise  ratio  based  on  the 
simplified  algorithms  of  Barton. 
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SECTION  I 


INTRODUCTION 

The  problem  treated  here  is  the  classical  detection  of  RF  target  signals  in  a Gaussian 

1 2 

noise  background.  This  was  initially  analyzed  and  presented  by  Marcum  and  Swerling  * 
and  extended,  since  that  time  by  a number  of  writers.  Since  tabulated  lists  and  curves  are 
often  somewhat  awkward  to  use,  the  writers  were  interested  in  hand-calculator  programs 
which  could  give  numerical  results  over  a wide  variety  of  detection  parameters.  In 
particular,  programs  for  HP-65  and  HP-67  calculators  are  presented,  but  the  algorithmic 
approaches  could  easil"  be  programmed  on  other  calculators. 

3 

Tne  writer  sta*.  this  work  by  following  Barton's  approach,  an  empirical  approxi- 

mation which  directly  provides  required  signal-to-noise  ratio  (SNR)  for  given  probability  of 
detection  (PD)  for  various  target  models.  The  accuracy  of  the  approximation  is  within  a 
dB  for  normal  parameter  values,  but  the  approach  does  not  lend  itself  to  finding  probability 
of  detection,  given  the  SNR,  without  an  iterative  approach  which  is  beyond  the  capability  of 
the  HP-65. 

4 

Later,  the  writer  found  an  excellent  report  by  Shnidman  who  had  found  finite  recur- 
sive series  solutions  for  probability  of  detection  for  the  basic  four  Swerling  target  models 
for  fixed-threshold  detection,  and  who  also  presented  infinite  series  algorithms  for  non- 
fluctuating and  the  generalized  chi-squared  target  models.  These  algorithms  weie  found  to 
be  directly  programmable  for  the  HP-65.  S.  P.  Applebaum  has  translated  these  programs 
for  the  HP-67  and  these  are  included  here.  For  the  special  case  of  a Swerling  Case  II 
target,  these  programs  can  solve  for  either  PD  or  for  SNR. 

In  spite  of  the  importance  of  both  CFAR  detection  and  noncoherent  integration,  there 
are  few  papers  in  the  literature  which  combine  the  two.  One  excellent  paper,  however,  is 
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that  of  Mitchell  and  Walker  . This  paper  treats  the  background  estimation  type  of  CFAR, 
and  it  has  provided  algorithms  used  here  to  cover  these  cases.  For  the  special  case  of 
Swerling  Case  II,  a finite  series  solution  and  programmed  iterative  inverse  can  provide 
either  PD  or  SNR  as  with  the  fixed- threshold  case.  For  other  target  models,  a truncated 
infinite  series  solution  provides  PD  given  the  SNR. 
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The  implementations  which  these  analyses  treat  are  illustrated  in  Figures  1-1  and  1-2. 
In  Figure  1-1,  the  noise  background  is  considered  to  be  constant  and  known.  N samples  of 
signal-plus-noise  are  summed  and  compared  to  a fixed  threshold.  If  the  threshold  is  ex- 
ceeded, a target  detection  is  declared.  The  threshold  is  set  according  to  a specified  false- 
alarm  probability  in  the  absence  of  signal.  In  the  recursive  solutions  programmed  here,  the 
first  step  is  the  calculation  of  the  threshold  value  given  the  false-alarm  probability. 

Figure  1-2  considers  the  case  where  the  background  noise  is  unknown  or  slowly  varying 
so  that  a noise  estimate  must  be  made  in  order  to  establish  a detection  threshold  value.  The 
noise  estimate  here  considered  is  the  sum  of  R independent  detected  noise  samples,  and 
corresponds  to  the  most  common  type  of  constant  false-alarm  rate  (CFAR)  detector. 

B6766 


Figure  1-1.  Fixed-Threshold  Detection 


RF  SIGNAL 
+ 

GAUSSIAN 
NOISE 1 


SUMMATION 

l 

SQUARE- LAW 

* 

OF  N 

V 

THRESHOLD 

DETECTOR 

SIGNAL ■(  NOISE 
SAMPLES 

1 'W 

DETECTOR 

DETECTION 

INDICATION 

fr. 


SUMMATION 
OF  R NOISE 
SAMPLES 


ADAPTIVE 

THRESHOLD 


-r  T 
(xa  ) 


Figure  1-2.  CFAR  Detection  (Background  Normalized) 
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The  target  models  used  here  follow  that  of  Swerling  and  the  integrated  target  SNR 
distributions  all  may  be  considered  special  cases  of  the  chi-squared  (or  gamma)  distribution 
given  by 

K _Kz 

W(z-Z)  - (5mTT  (I)  zK_1  6 Z <!-l> 

where  z is  the  integrated  target  power  SNR  for  any  trial  and  Z = z.  K is  a distribution 
parameter  which  can  have  any  value  greater  than  zero  and  certain  values  of  K correspond 
to  standard  target  fluctuation  cases*  as  follows: 

0 < K < 1 : Weinstock  case 

K = 1 : Swerling  case  I 

_ z 

w(z,  Z)  = | e Z (1-2) 

Exponential  power  or  Rayleigh  voltage  distribution.  Target  constant  over  N 
integrated  samples.  Results  from  radar  target  composed  of  many  separate 
scatterers.  Often  satisfied  by  aircraft. 

K = 2 ; Swerling  case  ill 
- 2z 

w(z,z)  = -^o  e Z 

Approximate  distribution  of  target  with  Rayleigh  and  fixed  components  of  equal 
average  power.  Target  constant  over  N integrated  samples. 

K = N : Swerling  case  n 

Same  basic  target  distribution  as  for  K = 1 but  with  target  amplitude  independent 
over  N integrated  samples.  Often  satisfied  by  aircraft  targets  with  pulse-to-pulse 
radar  frequency  agility. 

K = 2N  : Swerling  case  IV 

Same  basic  target  distribution  as  for  K = 2 but  with  target  amplitude  independent 
over  N integrated  samples. 

g 

* Nathanson  provides  ? good  discussion  of  these  target  models. 
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K = 00  : Nonfluctuating  target  or  case  0 


z = Z 

Note  that  other  values  of  K may  be  useful  for  cases  with  different  target  fluctuation 
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rates  or  with  block  correlation  within  the  N samples  integrated.  If  the  diversity 
order  within  the  N samples  is  Ne  then  with  a Rayleigh  target  K should  be  taken 
equal  to  Ne.  For  a Case  HI  target  distribution  K should  be  taken  as  2Ne. 


SECTION  II 


FIXED  THRESHOLD  DETECTION  - RECURSIVE  SOLUTION 


1.  SWERLING  CASE  II 

Although,  as  shown  later,  generalized  programs  can  be  written  to  cover  chi-squared 
distributions  of  any  K,  it  is  worthwhile  to  consider  some  special  cases  since  they  can  provide 
both  simpler  and  faster  running  calculator  programs.  We  shall  start  with  the  simplest  of 
these  cases  - fixed  threshold  detection  with  a Swerling  case  II  target  model. 

The  probability  density  function  of  v,  the  integrated  signal-plus-noise  variate,  is 
given  by,  * 


N-l 


f(v)  = 


V 

1+X 


>N 


(1+X)  (N-l)  ! 

where  X is  the  average  SNR  of  each  sample. 


(2-1) 

The  probability  of  detection  is  then  given  by 


P2 


■ / 


N-l 


_v 

1+X 


(1+X)N  (N-l)  ! 


dv 


(2-2) 


This  may  be  integrated  by  parts  to  give 
N-l 

P2  = 


m 


e 


Y 

1+X 


(2-3) 


— L m!  (1+X) 
m = 0 ' ' 

Notice  that  for  X = 0 this  reduces  to  the  false  alarm  probability  (PF) 
N-l 

PF  = 

m =0 


/ m! 


Y™  -Y 
e 


(2-4) 


~ Eqn.  IH.  10  of  Swerling  or  Eqn.  (39)  for  fN(V  |X)  of  Mitchell  and  Walker  on  page  675  noting 
that  our  X = Z/N  corresponds  to  their  X/N. 
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which  provides  an  implicit  solution  for  Y given  PF.  Notice  also  that  by  substituting 
Y - Y/l+X  in  Equation  (2-4),  we  obtain  Equation  (2-3)  so  that  one  program  routine  can  be 
used  for  both.  This  common  equation  shall  be  written  as 


m = 0 

with  each  definition  of  Y giving  the  appropriate  corresponding  definition  of  P. 

Gi/en  a desired  PF,  the  first  step  in  finding  either  P2  given  X or  X given  P2  must  be 
to  find  Y using  Equation  (2-5)  with  P = PF.  Perhaps  one's  first  thought  might  be  to  use 
Newton's  method  to  find  the  root  of  P - PIN  where  PIN  is  the  specified  value  and  P is  obtained 
from  the  equation  for  a given  Y.  However,  since  P versus  Y has  an  inflection  point,  con- 
vergence is  not  assured  and  it  is  better  to  use  In  (P/PIN).  This  leads  to  incrementing  Y for 
successive  trials  by 


AY  = - 


In  P-ln  PIN 


-pr  In  (P/ PIN) 


We  find  by  differentiating  Equation  (2-5)  that 


P'  = - 


(N-l) ! 


= - last  term  of  P series. 

til 

In  the  algorithmic  expressions,  the  m term  of  this  series  is  used  and  shall  be 
designated  YM.  Each  term  is  determined  recursively  from  the  previous  term.  After  com- 
pleting the  series,  we  will  have  in  storage  the  last  term,  YM,  so  we  can  use  this  for  -P'. 
Therefore,  in  applying  Newton's  method,  the  Y-increment  for  successive  trials  is  given  by 


P . P 

AY  YM  1 PIN 


(2-6) 


To  begin  the  iteration  of  Newton's  method  a starting  value  for  Y is  also  needed  which 
shall  be  designated  YO.  The  writer  found  empirically  that  the  .following  expression  approxi- 
mated Y quite  closely  for  small  values  of  PIN  as  may  customarily  be  desired  for  PF,  and 
was  programmable  with  very  few  program  steps  on  the  HP-65. 

YO  = N - VN  +2.3  sTl  ('/I!  + 's/N  - 1) 

L = - log  PIN  (2-7) 

Using  this  start  only  three  or  four  iterations  are  needed  to  calculate  Y to  10  significant 
figures  for  any  value  of  PIN  of  interest. 

This  solution  of  Equation  (2-5)  for  Y can  more  completely  be  specified  by  using  the 
algorithmic  notation  of  Iverson  (following  Shnidman's  practice)  as  given  in  Figure  2-1. 

A brief  explanation  of  this  notation  is  first  in  order.  The  arrow  notation  implies  a 
specification,  that  is,  the  statement,  L - — log  PiN,  is  translated  to  mean  that  the  quantity  L 
is  specified  by  - log  PIN.  The  normal  execution  of  the  statements  is  line  by  line  starting  at 
the  top,  but  a branch  may  be  designated  by  an  arrow  between  two  statement  lines.  A condi- 
tional branch  is  denoted  by  a colon  statement,  and  the  branch  is  executed  if  the  comparison 
condition  specified  oi.  the  arrow  is  safJsfied.  Otherwise  the  next  statement  in  the  sequence 
is  executed. 

The  brackets  labeled  D and  E on  Figure  2-1  correspond  to  subroutines  in  the  HP-65 

program  which  follows  and  are  shown  here  for  convenience.  Notice  that  the  iteration  is 
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terminated  when  I AY/Y|  is  less  than  10  . Since  the  stored  value  of  Y has  already  been 

corrected  by  the  indicated  AY  and  the  convergence  is  quite  rapid,  Y is  usually  accurate  to 
9 or  10  significant  figures. 

Having  obtained  Y for  a given  PF  using  this  algorithm,  we  can  calculate  P2  directly 
from  Equation  (2-5)  (i.e. , subroutine  E]  by  the  substitution  Y = Y/(l+X).  Alternatively, 
if  P2  is  given  and  it  is  desired  to  find  X,  we  can  substitute  PIN  = P2,  find  a corresponding 
Y2  using  the  program  of  Figure  2-1  an  then  find  X = Y/Y2- 1. 

These  features  are  all  contained  in  the  Program  HP-65  Y-P2  given  here.  Most  of  the 
program  comes  directly  from  the  algorithmic  program  of  Figuie  2-1,  but  a few  cc  uts 
may  help  in  its  understanding.  First  of  all,  the  writer  has  often  recordeo  in  the  v o ont 
space  on  each  line  of  the  HP-65  programs  the  stack  contents  in  the  order,  x,  y,  z,  t.  This 
may  be  useful  to  understanding  since  the  stack  is  often  used  in  these  programs  for  temporary 
storage.  This  practice  saves  on  use  of  the  storage  registers  which  is  sometimes  necessary 
and  also  often  leads  to  shorter,  faster  running  programs. 


YO  - N - /TP+  2.3  /T  (/T+  /N1-  1) 


Y «-  YO 
/ 

( 

^ YM  «-  e"Y 

YMS  «-  YM 

M <-  0 

M <-  M+l 

— N : M 

YM  «-  YM  x Y/M 
YMS  <-  YMS  + YM 

■>  P ■*-  YMS 


> D 


> E 


J 


AY  «•  (P/YM)  In  (P/PIN) 
Y «■  Y + AY 

•V  C 

10"°  : j AY/Y I 

I *-  EXIT 


j Figure  2-1.  Algorithmic  Program  for  Y Given  P 
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Next,  the  LBLE  subroutine  incorporates  LBLD  as  a loop  to  save  a program  step. 

When  the  D statement  is  reached  on  the  last  step,  it  jumps  to  LBLJ.  When  the  test  N < M 
is  satisfied,  the  RTN  jumps  back  to  the  E-subroutine  call,  since  the  HP-65  has  only  one 
program  step  register,  and  that  holds  only  the  initial  subroutine  call  step  number  no  matter 
how  many  successive  subroutines  are  called  before  the  RTN  statement  is  executed. 

Since  for  finding  X given  PF  and  P2,  the  Y-algorithm  programmed  here  is  used  twice, 
it  is  preceded  by  storing  the  previously  calculated  Y from  R2  in  R6.  The  first  time  through 
nothing  exists  in  R2  anyway,  so  these  steps  can  be  ignored  in  trying  to  understand  the  pro- 
gram. After  running  the  second  time  with  PIN  — P2  the  program  stops  with  Y2  displayed. 
Depressing  the  R/S  key  then  restores  Y in  R2  and  calculates  X from  the  Y2  and  Y previously 
found. 

Finding  P2  given  X uses  the  LBLA  function  and  is  straightforward,  requiring  only  that 
Y has  been  previously  calculated  or  otherwise  stored  in  R2. 

It  is  interesting  to  note  that  Shnidman  was  concerned  about  accuracy  of  the  calculation 

-Y 

and  underflow  for  certain  cases  such  as  e for  large  values  of  Y.  His  computer  was  equi- 
valent to  about  7 digit  words  and  he  went  to  double-precision  arithmetic  and  logarithmic 

±99 

calculation  in  underflow  cases.  With  the  10  digit  words  and  10  range  of  the  HP  calcula- 
tors, together  with  direct  monitoring  by  the  operator,  such  measures  are  really  unnecessary. 
The  programs  have  been  written  such  that  if  input  parameters  which  would  lead  to  underflow 
are  entered,  the  underflow  condition  results  almost  immediately.  This  .s  indicated  on  the 
HP-65  by  interruption  of  the  program  sequence  with  the  display  reading  zero. 

2.  SWERLING  CASES  I,  IH,  IV 

In  this  section,  the  basic  expressions  to  be  programmed  will  not  be  derived  but  will  be 
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taken  directly  from  Shnidman  to  which  the  reader  is  referred  for  more  detail. 

We  shall  refer  to  the  function  represented  by  Equation  (2-5)  as  P (N,  Y).  The 
probability'  of  detection  of  a Swerling  Case  I target  can  then  be  found  as 

PI  = e"Y^Z+1)  for  N = 1, 

or 


PI  = P (N-l,  Y)  + e~Y(Z+1)  £l  - p(n-1,  for  N>  2 

Here  the  integrated  power  signal-to-noise  ratio,  Z = NX  was  used. 


(2-8) 


2-5 


The  HP-6  5 PI  program  directly  implements  this  expression.  Prior  to  running  this, 
the  Y-P2  program  must  be  run  to  store  Y in  the  R2  register.  The  same  basic  LBL  E- 
subrcutine  is  used  here  as  in  the  prior  program,  modiPud  slightly  to  give  P(N-1,  Yl)  where 
Y1  is  in  the  X-register  prior  to  calling  the  subroutine.  The  coefficient  of  [ 1 - P]  is  also 
tested  so  that  if  too  small,  the  LBLE  subroutine  is  not  run  a second  time. 

In  a similar  way,  the  PD  for  a Swerling  Case  HI  target  can  be  found  as 


P3  = e 


(Z+2) 


for  N = 1 


P3  = 


YN~2  e~Y 
(N-2)  ! 


+ P(N-1,  Y) 


m 


N-2  - 


Mz2)  + _2Y‘|  . r _ 

Z Z+2 J [ 


P(N-1,  for  N>  2 


(2-9) 


The  HP-65  P3  program  directly  implements  this  expression.  It  is  run  following  the 
Y-P2  program  to  find  Y as  was  the  PI  program. 

Shnidman  shows  that  the  PD  for  a Swerling  Case  IV  target  can  be  written  as 

N - 1 „ 2N- 1 ,,  f M-N  v " 

P4  - y -V  r -V  1_'f'  N!  (_x\  /-M 

V ~ La  M!  La  m L^  K!(N-K)!  \X+2/  VX+2/ 


M =0 


2N- 1 

E 

e~V 

r M-N 

i-  E 

N! 

: / 2 \N  *K 

M! 

K ! (N-K) ! 

\X+2/ 

VX+2  / 

M=N 

K = 0 

J 

(2-10) 


where 


V = 2Y/(X+2). 

Although  the  first  term  is  clearly  equal  to  P(N,  V),  the  second  term  is  an  extended  summa- 
tion of  the  same  form  as  the  first  with  a more  complex  term-by-term  multiplier.  The  pro- 
gramming involves  a doubly  recursive  approach  which  is  best  illustrated  in  algorithmic 
form  on  Figure  2-2, 


i 

i 
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M 0 

ZK  - [2/(X+2)'JN 
V - 2Y/(X+2) 

YM  - e~V 
YMS  - YM 
M M+l  

— N : M 

YM  <-  YM  V/M 
YMS  - YMS  + YM  ■ 

SUM  - YMS 
ZKS  «-  ZK 


j 

| Figure  2-2.  Algorithmic  Program  for  P4 


4 

{ 

■f 
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For  relating  Figure  2-2  to  Equation  (2-10),  note  that  K = M-N  so  that  no  separate 
index  is  needed.  The  nomenclature  YM  corresponds  to  (VM/M ! ) e V while  ZK  corresponds 
to  each  term  of  the  K summation  and  YMS  and  ZKS  have  corresponding  relationships  to  YM 
and  ZK,  respectively. 

The  order  of  some  of  the  steps  listed  here  is  arbitrary  and  are  written  to  correspond 
to  the  program  HP-65  P4  for  consistency.  This  program  is  used,  as  for  PI  and  P3  after 
running  Y-P2  to  find  Y. 

3.  GENERALIZED  AND  NONFLUCTUATING  TARGET  MODEL 

For  the  nonfluctuating  target,  as  well  as  the  general  case,  the  summation  of  an 
infinite  series  is  required,  and  the  nonfluctuating  target  can  be  considered  a special  case 
of  the  general  formulation,  Mitchell  and  Walker  give  a straightforward  derivation  which  in 
our  nomenclature  can  be  written  as  follows . 

The  distribution  of  v for  a given  integrated  signal-to-noise  ratio,  z,  is  given  by 


fN(v|z)=(z)  e’(V+Z)lN-l(2'^) 


^ > yN+b-l  e-  (v+z) 
= b!  (N+b-1)  ! 


(2-11) 


as  given  by  Marcum  and  Swerling. 

The  probability  that  v will  exceed  a fixed  threshold,  Y,  is  then 


P(v  > Y |z.)  = j fN(v|z)dv 


N+b-1 


Ezr  -z  x_ 

b!  6 it /T 


(2-12) 


m =0 


For  a nonfluctuating  target,  this  gives  the  desired  PD  by  letting  z = Z.  For  a fluctuating 
target,  we  must  integrate  over  the  distribution  of  z as  follows 


00 

P = J w(z,  Z)  P (v  > Y|  z)  dz 


(2-13) 


Using  w(z,  Z)  for  the  generalized  chi-squared  distribution  of  Equation  (1-1),  this  yields 
the  generalized  PD, 


, N+b-1  „ 
b . — » „,m 


™ (K-fb-1)!  / _JK\  /Z\  V IZ  e"Y 

L^j  b ! (K-l)  ! \K+Z  / \K+Z/  Z-r  m! 


(2-14) 


m =0 


Shnidman  changes  the  order  of  summation  of  these  expressions  so  as  to  get  a more 
direct  measure  of  error  which  can  be  used  to  truncate  the  summation  to  a finite  number  of 


terms.  This  expression  then  becomes 


?G  = 


YM  + 


YM  1 


m =0 


M =N 


M-N 


-£ 


where 


YM  -Y 
YM  = M!  6 


XB  - 1 n_v)K  y = Z - 

A b!  (K-l) ! ' v'  v ’ K+Z 


(2-15) 


Note  that  a nonfluctuating  target  corresponds  to  the  limit  as  K — 00  for  which 

zb  -Z 
xb  - e z 


(2-16) 


The  error  in  Equation  (2-15)  for  a truncated  summation  is  shown  by  Shnidman  to  be 
given  by  the  product, 


m - 0 


M-N 


YM  1 


-£ 


(2-17) 
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The  programs  for  PG  given  here  test  this  product  after  each  term  of  the  PG  summation  and 
when  it  becomes  less  than  10  * , which  seems  a suitably  small  number,  the  summation  is 
stopped. 

An  algorithmic  program  for  PG  is  given  in  Figure  2-3,  the  program  corresponding 
directly  to  the  HP-65  PG  program.  As  with  the  other  programs,  it  is  necessary  to  run 
Y-P2  first  to  find  Y.  The  initiation  of  the  PG  program  requires  entry  of  both  X,  and  the 
target  distribution  parameter,  K. 

Suitable  values  of  K were  discussed  in  the  introduction,  but  the  special  case  of  a non- 
fluctuating target  provides  some  difficulty  since  infinity  is  not  an  allowable  entry  value. 

5 

The  best  large  number  to  substitute  for  infinity  in  this  case  was  found  to  be  about  10 
(entered  with  only  two  keystrobes  as  EEX  5).  Larger  values  for  K give  difficulties  for 
3ome  values  of  Z = NX  in  calculating  the  initial  value  of  XB,  while  smaller  values  are  less 
accurate  approximations  of  infinity.  This  compromise,  however,  apparently  gives  ar*  ac- 
curacy for  the  calculated  P0  of  at  least  three  places  for  any  value  of  Z.  To  avoid  the 
required  entry  of  K in  this  case,  as  well  as  to  provide  greater  accuracy  if  wanted,  a 
modified  form  of  the  PG  program  is  given  here  as  the  P0  program  which  calculates  the 
detection  probability  for  a nonfluctuating  target.  This  is  based  on  using  Equation  (2-16) 
for  XB  in  place  of  the  more  general  Equation  (2-15).  This  eliminates  the  computation 
difficulty  for  a large  value  of  K,  since  it  is  not  used. 
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_Y 

YM  - e T 
YMS  «-  YM 
M «-  1 


N : M 

YM  - YM-Y/M 
YMS  «-  YMS  + YM 
M «-  M+l  

SUM  - YMS 
V <-  Z/(K+Z) 
X8  «-  (1-V)K 
ZK1  «-  K 
XBS  <-  XB 


A 


f 

t 


YM  - YM  Y/M 

YMS  «•  YMS  + YM  l,  A 

SUM  <-  SUM  + YM  (1-;:bs)  f 

M *■  M+l  ) 

XB  «-  XBVZK1/(M-N) 

ZK1  «•  ZK1+1 

(■  : ( 1 -VMS ) fl-XBS) 

p SUM  ► EXIT 


Figure  2-3.  Algorithmic  Program  for  PG 
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SECTION  III 


FIXED  THRESHOLD  DETECTION  - BARTON  ALGORITHM 
3 8 

Barton  and  Cann  were  interested  in  a somewhat  universal  set  of  curves  which  could 
be  used  simply  to  find  radar  detection  performance  over  various  target  and  radar  parameters. 

They  found  that  an  ideal  detector  curve  plus  a set  of  relatively  simple  Ids*  factors,  i.  e. , 

detector  loss,  integration  loss,  coUapsing  loss,  and  fluctuation  loss,  gave  very  reasonable 

accuracy  for  normal  values  of  PD,  PF,  and  N,  and  for  the  target  distributions  we  have  been  ; 

considering.  Barton's  algorithms  have  also  been  programmed  for  the  HP-65  and  HP-67  and 

are  included  here.  These  programs  are  complementary  to  the  PD  programs  previously 

given  in  that  they  calculate  a required  signal -to-noise  ratio  for  a given  PD  rather  than  the 

other  way  around. 

For  the  nonfluctuating  target,  and  no  collapsing  ratio,  the  Barton/'Cann  algorithm  can 
be  written 

SNR  (dB)  - 10  log  j | [X  + n/X(X+ 9.2) 
r 'i  2 

x = ^[Q'VfJ+Q-VdiJ 
°°  _ tf. 

Q(y)  = f e 2 dt  (3-1) 

247  ' 

(To  correlate  this  with  Barton's  nomenclature  Xq  = 2X) 
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The  inverse  Q function  is  calculated  by  the  approximation.  For 
. . a + a.t 

P<  - , Q“  (P)  = t - — o (3"2) 

z l+b^+bgt 


3-1 


where 


ax  = 0.271, 
bx  » 0.992, 
b2  = 0.0443, 

t = V In  (l/P2)  , 

and  for  P > 1/2,  Q-1(P)  = -Q“1(l  - P)  . 
Also,  integration  gain  is  given  by: 
Gi  ==  SNRN  - SNR1 


These  are  directly  programmed  in  the  HP-65  SNRN  and  HP-67  SNR  programs. 

For  a nonfluctuating  target  with  collapsing  loss,  or  for  a fluctuating  target,  the  HP-65 
SNRF  and  HP-67  SNR  programs  make  the  following  calculation: 

SNRF  (dB)  = SNRN  + log  —■  + Lf  , (3-3) 


where  Lf,  the  fluctuation  loss,  is  given  by 

um  . lojsKSjLmk 


(3-4) 


SNRN  is  the  value  calculated  by  program  SNRN  for  N and  SNR1  is  the  value  calculated 
by  program  SNRN  for  N = 1.  D is  the  single  pulse  average  SNR  for  fluctuating  target  detec- 
tion and  depends  on  the  target  model  used.  For  the  Rayleigh  fluctuation  model  of  Swerling's 
Cases  1 and  2 


In  PF  _ 
In  PD  “ 


(3-5) 


i 


i 


i 

» 

t 

t 

\ 

i 


For  the  one  dominant  plus  Rayleigh  fluctuation  model  of  Swerling’s  Cases  III  and  IV,  D . is 
10 

given  implicitly  by 


PD  = 


2D34lnPF  j pF-  “34 


<2  + D34> 

The  writer  found  that  the  solution  to  this  equation  is  well  approximated  by 
D„.  = (0.361  - log  PD)  ( — — 7 - - 1.29  - 0.96  VI  - PD  1-2 

34  \ sTTpb  I 


(3-6) 


(3-7) 


and  this  expression  is  used  in  HP-35  SNRF  and  H >-67  SNR  to  avoid  the  need  for  reiteration. 
The  greatest  error  in  this  approximation  occurs  for  low  values  of  PD  but  it  is  accurate  to 
better  than  0. 5 dB  ior  PD  equal  50%  and  within  1 dB  for  PD  equal  30%.  1'or  PD  greater  than 
90%,  it  is  accurate  to  within  0. 2 dB.  An  extremely  bad  choice  of  PD  too  low,  or  PF  too  large 
may  cause  D34  to  be  negative  and  flashing  zeros  will  indicate  this  error  when  running  the 
program. 

Finally,  the  programs  calculate  the  diversity  gain  as 
Gd  (dB)  = (N  - 1)  Lf 

V/ 

and  the  range  ratio,  re  Swerling 
SNRF/40 

R/R  = 10 
o 


f 


J 
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SECTION  IV 


CFAR  DETECTION  - RECURSIVE  SOLUTION 
1.  SWERLING  CASE  II 


As  with  the  fixed  threshold,  the  Case  II  target  model  leads  to  a simple  analysis  and 
finite  summation  for  finding  the  probability-of-detection.  Starting  with  Equation  (2-3),  the 
constant  Y can  be  replaced  by  the  variable  y to  have 

m 


N-l  ( __ JL] 
l x+1  j 

p<v>y>  = L mi 

m = 0 


_ JL_ 
X+1 


(4-1) 


for  each  specific  value  of  y. 


Note  on  Figure  1-2  that  y is  derived  from  u and  that  u is  the  sum  of  R independent 
Rayleigh  noise  samples  of  unit  average  power  - unity  since  we  also  normalized  the  magnitude 
of  v to  the  average  noise  power.  Therefore,  u has  the  distribution 


R-l 


p(u)  = 


-u 


(R-l)! 

Then  the  overall  probability  of  v exceeding  y is  given  by 


(4-2) 


P = j p(u)  P (v  > ^ j du 
0 


(4-3) 


where  T is  a calibrating  factor  which  must  be  set  to  achieve  the  desired  false -alarm 
probability  and  is  analogous  in  our  further  derivation  here  to  Y which  determined  the  faise- 
alarm  probability  in  the  fixed  threshold  case.  Substituting  Equations  (4-1)  and  (4-2)  into 
Equation  (4-3),  interchanging  the  order  of  summation  and  integration,  and  integrating, 
one  gets 


N-l 

N-l 

D 

E ^ 

m = 0 

•E 

m=0 

(R+m-1)  ! 

(T2) 

(4-4) 

m!  (R-l)! 

(T2+l)R+m 

where  T2  = T(X+1). 


/ 
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In  a similar  manner  to  that  for  finding  Y previously,  let  X = 0 so  that  T2  = T and  find 
the  varne  of  T for  which  P equals  the  false-alarm  probability. 

This  process  is  best  done  by  using  Newton's  method  on  In  (P/ PIN)  as  before  so  that 

AT  = (4-5) 

dT 


and  we  find  from  Equation  (4-4)  that 


dp 

dT 


R 

T 


m = 0 


E 


m P 


' m 


(4-6) 


Denoting  the  last  summation  as  Q,  Equations  (4-4),  (4-5)  and  (4-6)  yield 

„ lnip/nm 

q/p  - R/T  1 

Since  the  terms  of  Q are  closely  related  to  those  of  P,  both  sums  can  be  formed  at  the 
same  time. 

We  are  left  with  the  problem  of  the  initial  value  to  use  for  T.  Extending  the  curve 
fitting  approach  of  before  the  writer  found  a reasonable  initial  value  to  be  given  by 


_1_ 

TO 


1 - (KN)1/R 


(1-B)  N+B 


/ N -n/n 

l 2.3L 


, ^L  + Vn-1  \ 
\/"l  ) 


R-l 

R+0. 922  ’ 


-log  PIN 


(4-3) 


The  value  of  TO  from  Equation  (4-8)  was  found  to  provide  a sufficiently  good  start  foi  itera- 
tive convergence  over  the  range  of  10  10  < P < 1 and  1 < R < 1000. 

Unfortunately,  this  takes  more  than  50  program  steps  so  that  a separate  program  card 
is  necessary  for  data  entry  and  calculation  of  TO  with  the  HP-65.  After  running  this  program, 
HP-65  P2C  can  be  used  to  perform  the  iterative  calculation  of  T for  a given  false-alarm 
probability  and  P2  for  various  input  values  of  SNR  = 10  log  X.  If  SNR  is  to  be  found  for  a 
given  P2  program,  HP-65  TO  must  be  rerun  with  P2  input,  followed  by  HP-65  P2C  again, 
the  process  being  directly  analogous  to  that  of  the  HP-65  Y-P2  program  for  a fixed  threshold. 
Since  the  HP-67  has  mere  program  storage,  the  TO  calculation  is  included  in  the  HP-67  P2C 
program. 
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2.  GENERALIZED  TARGET  MODEL 


The  relationship  of  the  CFAR  process  to  be  fixed  threshold  process  in  general  is  the 
same  as  it  was  for  Case  II.  Starting  with  Equation  (2-14),  let  Y = u/T  and  integrate  over  the 
distribution  of  u from  Equation  (4-2)  to  find  the  overall  PD.  By  this  process,  one  obtains 

co  N+b-1 

XB  ^ PM 
b=0  m=0 


An  algorithmic  program  for  PGC  is  given  in  Figure  4-1  which  follows  closely  that  for 
PG.  Unfortunately,  this  could  not  be  fitted  into  100  HF-65  program  steps  so  it  had  to  be 
programmed  on  two  cards,  HP-65  PGC(l)  and  PGC(2).  PGC(l)  incorporates  the  iteration 
for  T given  TO  and  finds  A = 1/(T+1).  Therefore,  PGC  is  found  by  running  in  sequence 
HP-65  TO  to  find  TO,  PGC(l)  to  find  A and  to  enter  SNR  and  K,  and  finally  by  running  PGC(2) 
as  noted  on  the  CFAR  Detection  instruction  sheet.  For  the  HP-67,  the  entire  calculation  is 
included  on  the  single  program  card,  HP-67  PGC. 

3.  CFAR  LOSS 

The  increase  in  SNR  required  with  a CFAR  detector,  as  compared  to  a fixed  threshold 
detector,  has  been  termed  CFAR  loss.  This  concept  is  a convenient  one  because  the  CFAR 
loss  is  essentially  independent  of  the  target  fluctuation  model,  at  least  as  far  as  the  five 
Marcum  and  Swerling  models  are  concerned.  Although  many  papers  in  the  literature  deal 

g 

with  CFAR  loss  for  various  CFAR  detector  schemes,  the  paper  by  Mitchell  and  Walker  is 
the  only  one  found  by  the  writer  to  cover  the  combination  of  noncoherent  signal  integration 
with  a background  normalizer  threshold.  Using  the  HP-65  Programs  of  this  paper  for 
Case  II  targets,  the  data  of  Figures  4-2  through  4-5  were  calculated  and  are  presented  here 
for  convenience.  The  loss  values  from  these  curves  may  be  used  as  a.  correction  to  fixed 
threshold  SNRs  for  the  fixed  threshold  programs.  This  may  be  handier  than  running  the 
CFAR  programs  for  many  cases  because  of  the  long  running  time  of  the  generalized  CFAR 
program.  Although  the  writer  has  verified  in  a few  sample  cases  that  other  target  models 
and  other  detection  probabilities  give  essentially  the  same  CFAR  loss  values,  it  will  be  left 
as  an  exercise  for  the  reader  to  be  convinced  that  this  is  true  for  the  cases  of  concern. 
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PM  <-  (1-A)R 
PMS  <-  pm 
M «-  1 


- N : 1 

PM  4-  PM.  B+JhI.  A 
M 

PMS  4-  PMS  + PM 
M 4-  M+l  

► V 4-  Z/(K+Z) 

XB  4-  (1-V)K 

XBS  4-  XB 
SUN  4-  PMS 


pm  4-  pm  BJhI  # A 

M 

PMS  - PMS  + SUM 
M 4-  M+l 


SUM  4-  SUM  + PM  ( 1 -XBS ) 
XB  ^ XB  -^N'1  .V 


XBS  4-  XBS  + XB 


i : ( 1 - XBS ) ( 1 - PMS ) 
p «.  SUM  


> 


A 


EXIT 


Figure  4-1.  Algorithmic  Program  for  PGC 
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B6769 


Figure  4-4.  Square-Law  CFAR  Loss  Swerling  Case  II  Target,  PF  = 10  PD  =0.5 


SECTION  V 


HP-65  PROGRAMS 

1.  FIXED-THRESHOLD,  RECURSIVE  SOLUTIONS 

These  HP-65  programs  calculate  the  probability  of  detection,  given  the  number  of 
samples  noncoherentiy  integrated,  the  false-alarm  probability,  and  the  average  sample  signal- 
to-noise  ratio,  for  the  various  Swerling  target  models.  The  PG  program  does  this  for  the 
generalized  chi-squared  target  model.  The  Y-P2  program  must  be  used  to  calculate  the 
threshold  value,  Y,  before  using  any  of  the  other  programs,  and  in  addition  it  can  calculate 
for  a Case  II  target,  either  probability  of  detection  given  average  signal-to-noise  ratio  or 
average  signal-to-noise  ratio  given  probability  of  detection. 


Specific  user  instructions  are  as  follows: 


STEP 

INSTRUCTIONS 

INPUT 

DATA/UNITS 

iveYS 

OUTPUT 

DATA/UNITS 

1 

Enter  program  Y-P2 

N 

|STO|  [T 

i 

PF 

|A  )f 

i 

Y 

Go  to  step  2,  3,  4 or  6 

1 II 

i 

2 

For  P2  (Repeat  or  go  to  step  3 as  desired) 

X 

IB  I L . 

! 

P2 

3 

For  X given  P2 

P2 

lA  . J f.  ' 

i 

Y2 

Ir/s  I i 

1 

X 

Repeat  or  go  to  step  2,  4,  or  6 as  desired 

f ‘ I 1 .... 

i 

_ j 

4 

For  PO,  PI,  P3,  or  P4  enter  that  program 

HE 

n 

5 

X 

Q| 

H 

D 

Repeat  or  go  to  step  4 or  enter  Y-P2  and  go  to 

i a 

H 

step  2 as  desired 

H ■ 

n 

6 

For  general  target  model  enter  program  PG 

i 1 1 ; 

1 

7 

K 

i i rr. 

1 

X 

IB  If'. 

i 

j 

P 

Repeat  or  go  to  step  4 or  enter  Y-P2  and  go  to 

r ] i : ■ 

I 

step  2 as  desired 

-l—l  l— 

i 

-j — 

i 


K 

N 

Ne 

PF 

P 

PO 

P1-P4 

X 

Y 

Y2 


Chi-squared  distribution  parameter  

Signal  and/or  noise  samples  integrated 

Target  diversity  within  N samples 

Probability  of  false  alarm 

Probability  of  detection  for  chi-squared  target 

Probability  of  detection  for  nonfluctuating  target 

Probability  of  detection  for  Cases  I-IV 

Avg  sample  power  S/N  within  N samples 

Fixed  detection  threshold 

Y/(l+X) 


f Case  0 : K=105  (~3  place  acc) 

Case  I : K=1 

Case  n : K=N 

(Genl  Rayleigh  target  : K=Ne) 

Case  HI  : K=2 
Case  IV  : K=2N 

(Genl  Rayleigh  + equiil  constant  target  : K=Ne) 
Weinstock  : 0 < K < 1 
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HP-65  Program  Form 

Titk  HP-65  Y-P2  Program  Listing Page of 


SWITCH  TO  W PRGM  PRESS  [f  ] PRGM  • TO  CLEAR  MEMORY 


j KEY 

1 ENTRY 

CODE 

SHOWN 

COMMENTS 



KEY 

ENTRY 

CODE 

SHOWN 

~ 

COMMENTS 

1 

REGISTERS 

ST04 

3304 

PIN 

1 

01 

R i N 

RCL2 

3402  '( 

- SAVE  y IN  R6 

RCL2 

3402 

Y 

ST06 

3306  / 

R/S 

84 

RCL4 

3404 

PIN 

RCL6 

3406 

RESTORE  Y IN  R2 

R,  Y 1 

7 

31 

ST02 

3302  7 

LOG 

08 

_E*nz1 

3507 

Y2,  Y 

CHS 

42 

L ~ ' * 

_«  . 

r3 

' RCL1  4 

3401 

N 

l 

01 

7 

41 

- 

51 

t 

31 

R/S 

84 

v - v / ^ 1 
x/  "■  1 

R*  PIN 

09 

n/N.  N.  L 

LBL 

. 23  _ 

i 

51 

N-n/n.  L 

B 

12 

X _ . 

i 

IgLSTX 

3500 

n/N 

1 

01 

R5 

1 

01 

61 

- 

51 

«/T-l,  N - ,7ni  L— ■ 

) RCL2  r 3402 

SAVE  Y IN  R6 

m. 

3509 

L 

ST06 

3306 

Re  Y i 

* 

31 

gx<?y 

3507 

TT 

09 

\Tl 

~ 

81 

Y2  -Y/(l+X) 

+ 

61 

'*TL+*rS-it  n -n/n 

ST02 

3302 

R7  ^ 

gLSTX 

3500 

r \[l 

E 

15 

P2 

X 

71 

RCL6 

~3406 

RESTORE  Y IN  R2 

2 

02 

ST02 

3302 

Rh  YMS 

• 

83 

gRI 

3508 

i 

3 

03 

R/S 

84 

P2 

X 

"2,3  VL(s/lW?W£  N-n/N  ' 

LBL 

23 

Rs  

+ 

61 

E 

15 

ST02 

3302 

Y 6 

! RCL2 

3402 

Y 

LBL 

23 

CHS 

42 

LABELS 

1 

01 

f-1 

32  1 

A 

£ 

15 

p 

LN 

07 

~YM—  e"Y 

R X-P2 

gRt 

3509 

YM 

ST08 

3308 

YMS  — YM 

C . J 

-r 

si 

P/YM 

0 

00 

M~-0 

D YM  Loop 

1 RCL8  | 3408 

P 

LBL 

23 

F Y-P 

RCL4 

3404  . HN 

.P  . 

14 

M,  YM 

0 

-r 

81 

P/HN,  P/YM  ] 

1 

01 

1 

f 

31 

*- 

61 

M-  M+l 

2 

LN 

07 

AY— (P/YM)  LN(P/HN) 

RCL1 

3401 

N 

3 

X 

71 

3522 

NsM 

4 

STO 

33  : _ 1 

RCL8 

3408 

P-YMS 

5 

■if 

61 

RTN 

24 

6 

2 

02 

(Y  - Y+AY) 

RCL2 

3402 

Y,  N,  M,  YM 

7 

RCL2 

3402 

Y 

gRt 

3509 

YM.  Y.  N.  M 

8 

-r 

81 

x 

_ 71. 

YM.  Y.  N.  M— 

9 

_g 

35 

gRt. 

3509  1 

M.  YM.*  Y,  N.  M 

ABS 

06 

lAY/Yl 

81 

YM-YM  • Y/M 

FLAGS 

EEX 

43 

Tsto 

33 

1 

CHS 

42 

+ 

._61 

6 

06 

-6 

8 

08 

(YMS  -YMS  + YM) 

2 

gxs^y. 

3522 

10  S lAY/Yl 

.gRL 

3509 

OTO 

22 

D 

14 

M,  YM 

-■ 

v T0  lW  PROGRAM  INSERf  MAGNETIC  CARD  WITH  SWITCH  SET  AT  W/PRGM 
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HFUrLt. 


HEW  L 


HP-65  Program  Form 

Tl(1  HP-65  PI  Program  Listing  


Pnge of 


SWiTCH  10  W PRGM  PRESS  f ' [ PRC*M  j TO  CLEAR  MEMORY 


KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

REGISTERS 

RCL1 

3401 

'n  L_  f x 

E 

15 

Ri  N 

X 

71 

\ B 

ST07 

3307 

Y1 

ST03 

3303 

z 

CHS 

42 

RCL2 

3403 

Y 

fl 

32 

R?  Y 

gx— y 

3507 

Z.  Y 

LN 

07 

YM  - e"Y1 

1 

01 

ST  08 

3308 

YMS  — YM 

+ 

61 

Z+l,  Y 

1 

oi 

M— 1 

r3  Z 

ST04 

3304 

LBL 

23 

T 

81 

D 

14 

IcHS 

42 

RCL1 

3401  1 

N 

r4  Z+l 

ri 

32 

1 

01 

Z7(Z+1) 

LN 

07 

- 

51 

N-l 

R, 

ST05 

3305 

_ . . 

e-Y/(Z+l) 

_gx^y 

3522 

N-l  £ M,  YM 

RCL1 

3401 

N 

RCL8 

3408 

P— YMS,  N-l,  M.  YM  1 

s 

1 

01 

RTN  J 

24 

gx=y  | 

3523 

1=N 

RCL7 

3407 

Yl,  N-lt  M,  YM  ^ 

BePlN-l.Y) 

RCL5 

3405 

gRt  ... 

3509 

YM 

R/S 

84 

PI— e~Y/(Z+l) 

X 

. 71 

YM'Yl 

RCL2 

3302 

Y 

_*R  t 

3509 

M 

R7  Yl 

,E 

15 

PIN-1.  Y).  N-l 

T 

Fl 

YM-  YM  • Yl/M 

ST06 

3306 

STO 

33 

RCL3  1 

3403 

Z 

+ 

61 

iRa  YMS 

'RCL4 

3404 

Z+l,  Z.  P,  N-l 

8 

08 

(YMS -YMS  + YM) 

+ i 81 

' gR  t 

3509 

M 

ST04 

3304 

Z/(Z+1) 

l _ 

01 

Rq  Used 

gRt 

3509 

N-l 

+ 

61 

M— M+l 

g 

35 

D 

14 

ZX 

05 

I Z/tZ+ltt1*"1 

LABELS 

RCL5 

3405 

A _ 

V 

81 

hi 

- - 

B 

[ ST05 

3305 

S-  [ Z/(Z+1)]  N-le- 

■T~  " 

C 

fEEX  , 43 

n YM  Loop 

7 

07 

p Y—P 

gxs  y 

3522 

107==S 

0 

RCL6 

3406 

L ..  _ . 

1 

[R/S  • 84 

PI—  P(N-1,  Y) 

2 

CLX 

44 

'*  

3 

PRCL4 

3404 

Z/lZ+l) 

4 

RCL2 

3402 

Y 

5 

X 

71 

Vz7(z+ij 

<r 

6 

E 

1 15 

Pl-N-l,  YZ/(Z+1)1 

7 

CHS 

42 

ft 

1 

01 

3 

+ 

61 

1-P 

RCL5 

3405 

s 

FLAGS 

-r 

61 

(1-P)/S 

i 

1 

RCLb 

340t 

PIN-1.  Y> 

+ 

61 

2 

r7s 

84 

P 

LBL 

23 

i ' 



TO  RF-  'RD  p (OOHAM  MUSE*-  T MAGNE  TIC  CARO  W'TM  SWITCH  SET  AT  W/PRGM 


HEWLETT  JUT  PACKARD 


\ 

i 


HP-65  Program  Form 

T»tie  HP-65  P3  Program  Listing Page ot 


SWITCH  TO  W PRGM  PRESS  f j [_PRGMr]  T0  CLEAR  MEMORY 


KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS  j 

REGISTERS 

RCL1 

3401 

N L <x 

RCL5 

3405 

2Y/(Z+2) 

Ri  N 

X 

71 

\ B 

RCL3 

3403 

Z/2 

2 

02 

X 

71 

YZ/(Z+2) 

— 

81 

E 

15 

P(N-1  , YZ7(Z+2)) 

R?  Y 

ST03 

3303 

Z7~2 

CHS 

42 

RCL3 

3403 

i 

01 

1 

01 

+ . 

61 

1-P 

R3  Z/2  H 

+ 

61 

l+Z/2,  Z/2 

RCL5 

3404 

C 

-r 

81 

Z/(Z+2) 

X 

71 

C(l-P) 

RCL2 

3402 

Y 

ST04 

3304 

r4  c “• 

klstx 

3500 

. l+Z/2*  Y . .. 

RCL2 

3402 

P 

C(l-P) 

■f 

81 

r 2 Y/(Z +2)  , Z/(  Z +2) 

E 

15  ~ 

PiN-i  Y)  _ 

ST05 

3305 

1 ER  i 3509 

YM  =YN"2e“Y/(N-2)  1 

R52Y/(Z+2) 

RCL1 

3401 

N 

RCL5 

3405 

2Y/(Z+2) 

2 

02 

X 

71 

3522 

2 £ N 

+ 

L_61_ 

P(Y)  + YM  • 2Y/(Z+2) 

R6  N-2 

GTO 

22 

RCL4 

3404 

C(l-P) 

1 

01 

+ 

61 

kR  t 

3509 

R/S 

84 

P3 

R7  Y3~ 

«Rt 

3509 

LBL 

28 

X 

71 

2YZ/(Z+2)JS 

E 

15 

1 

01 

ST07 

3307 

Y3 

R0  YMS 

+ . 

-61  .. 

CHS 

. ' 42  . 

RCL5 

3405 

2Y(Z+2) 

f-i 

32 

CHS 

42 

LN 

07 

YM  — e~Y3^ 

Rg 

if-1.  i __J2  _ 

ST08 

3308 

YMS— YM 

LN 

07 

L.  -2Y/(Z+2)  _ nl 

1 

01 

M— 1 

X 

71 

e [l+2YZ/(Z+2>] 

LBL 

23 

LABELS 

R/S 

84 

D 

14  J 

A 

LBL 

23 

Trcli 

3401 

N 

B 

1 

01 

1 

01 

C 

- 

51 

N-2 

- 

51 

N-l 

D YM  Loop 

ST06 

3306 

gx-y 

3522 

N-lsM,  YM 

F Y-P 

eR  t 

3599 

Z/(Z+2).  N-2,  2Y/(Z+2) 

RCL8 

34  08 

P—YMS.  N-l.  M.  YM 

n 

f 

31 

RTN 

24 

1 

LN 

07 

RCL7 

3407 

Y3,  N-l,  M,  YM 

2 

X 

71 

-(N-2)  E((Z+2)/Z] 

IgRt 

3509 

YM 

3 

+ 

61 

EXP-2 Y/(Z+2)-(N-2)lnf  J 

X 

"71" 

YM*  Y3 

4 

CHS 

42 

ER  t 

3509 

M 

5 

r1 

32 

81 

YM-  YM  • Y3/M 

6 

LN 

07 

e-EXP 

STO 

33 

7 

1 

01 

+ 

61 

8 

RCL6 

3406 

8 .... 

08 

( YMS — YMS+YM) 

9 

RCL3 

3403 

Z/2.  N-2. 1.  e_EXP 

J?R  t 

3509 

M 

• 

81 

2(N-2)/Z 

i H 

01 

FUGS 

- 

51 

l-2(N-2)/Z 

+ 

61 

M-M+l 

1 

RCL5 

3405 

2Y/(Z+2)  ...  _ 

D 

14 

+ 

61 

2 

x 

71 

C— [l-2(N-2)/z+2Y/  (Z+2)]  t 

"LAr 

T 

ST  04 

3304 

1 

lioo 

q3?P-06l6  TO  RECORD  PROGRAM  INSERT  MAGNE-  TIC  C ARC  WITH  SWITCH  SET  AT  W/PRGM 


HP-65  Program  Form 

T tic  HP-6S  P4  Program  Listing 

SWITCH  TO  W'PRGM  PRESS  T)  f PRGM  TO  CLEAR  MEMORY 


Page of 


HP-65  Program  Form 

Tale  HP-65  PG  Program  listing 


Page of 


SWITCH  TO  W PRGM  PRESS  f ' [ PRGM^  10  ClCilR  MEMORY 


KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

REGISTERS 

RCL1 

3401 

n,x;k.  f 1 K 

RCL3 

3403 

V 

R,  N 

X 

71 

. _ ru-.i  t 

X 

71 

XB  * V/B 

ST03 

3303 

2 

T x 

gx— y 

3507 

ZK1 

gRi  1 

3508 

* B 

X 

71 

R,  Y 

ST04 

3303 

K 

ST04 

3304 

XB-XB  • ZK1  • V/(M-N) 

RCL2 

3402 

1 

gLSTX] 

3500 

ZK1 

CHS 

42 

IgRI  , 

3508  1 

XB,  (1-XBS),  (1-XBS),  ZK1 

r3  ^ 

r1 

32 

- 

51 

(1-XBS)— (1-XBS)  - XB 

V 

LN 

07 

gRt 

3509 

ZK1.  (1-XBS) 

ST07 

3307 

YM-e’ 1 

1 

01 

R4  k 

ST08 

3308 

YMS— YM 

+ 

61 

ZK1— ZKl+1 

XB 

1 

01 

■f 

L 

oi 

ST06 

3306 

M— 1 

j 

RCLSI 

3408  j 

YMS,  1,  ZK1,  (1-XBS) 

_ SUM  ~ 

LBL 

23 

1 

i - 

51 

1 

01 

1 

Vt 

3509 

1 

RCL1 

3401 

N,  M 

X 

71 

(l-XBSWl-YMS).  ZKl.d-J 

¥ 

s 

L 

3522 

N £ M 

EEX 

43 

GTO 

22 

CHS 

42 

2 

02 

( YM-YM-  Y/M 

8 

08 

io~8.  on 

Ry  YM 

A 

.11 

M • YMS- YMS  + 1M 

lfflc>y  . 

.3424 

GTO 

22 

( M-M+l 

1 

[RCL5 

3405 

rP— SUM 

'i 

01 

i 

R/S 

84 

R8  YMS 

LBL 

23 

t 

3509 

2 

02 

gRt 

3509 

ZK1,  (1-XBS) 

RCL8 

3408 

YMS 

GTO 

22 

Ko  logic 

fST05 

3305 

SUM— YMS 

3 

03 

i'“ 

' 01 

LBL 

23 

RCL3 

3403 

Z 

A __ 

h lL 

f” ZK1,  (1-XBS) 

LABELS 

RCL3 

3403 

z 

RCL7 

34)7 

YM 

a Recursloii 

RCL4 

3404 

K,  Z,  Z,  1 

RCL2 

3402 

Y,  YM,  ZK1,  (1-XBS) 

B . 

1 + 

61 

X 

71 

YM*  Y 

C . 

It 

[ si  u 

__ 

RCL6 

3406 

M.YM-Y,  ZKl.a-:iBS) 

D _ 

ST03 

3303 

r V— Z/(K+Z) 

T 

81 

E _ 

- 

51 

1 

-V 

i 

ST07 

33  17 

YM-YM*  Y/M 

O 

RCL4 

3404 

r k 

STO 

i3 

1 Used 

g 

35 

L 

,±.  .... 

61... 

? Used 

05 

r 

8 

_ . 08 

(YMS- YMS  + YML 

3 Used 

ST04 

3304 

xB-ii-vr 

X.Rt 

3509 

fl-XBS).  YM.  ZKlll-XBS) 

4 

gLSTx 

3500 

r ZKl—K.XB,  1— 

X 

71 

, YI4(l-XBS)xZKl,X)  - 

5 

IgRI 

3508 

XBS-XB,  1, 1,  ZK1 

STO 

33 

6 

51 

G-XBS),  1,  ZK1- 

+ 

61 

7 

gRI 

3509 

ZK1,  (1-XBS) 

5 

. _«L 

(SUM— SUM+YM(1-XBS)) 

8 

LBL 

23 

J 

YM— YMS*  Y/M 

! 

jjRl 

3505 

LZ  K1,  n . 1 

9 

3 

l~  03 

i YMS— YMS+YM 

1 

f 

RCL6 

3400 

M 

A 

11 

M1 

SUM-SUM+ YM(  1-XBS) ! 

1.  . 

01 

1*M,  ZKl^l-XBSl 

FLAGS 

RCL1 

3401 

N 

M-M+l 

+ 

61 

1 

- 

51 

M-N 

ST06 

3306 

M-M+l 

RCL4 

3404 

XB 

RTN 

24 

2 

gx — V 

3507 

m-N.XB.  ZK1.  (1-XBS) 

4} 

31 

XB/(M-N) 

V 'O 

5-7/5-8 


f4»K> 


TO  RECORD  PROCRAM  INS-RT  MAGNETJC  CARD  WITH  SWITCH  SET  AT  W PROM 


2. 


FIXED-THRESHOLD,  BARTON  ALGORITHM 


These  HP-65  programs  calculate  the  required  average  SNR  in  dB  for  a given  detection 
probability  and  target  model.  The  target  model  is  specified  by  its  probability  density  function, 
i.e.,  nonfluctuating,  Rayleigh  (ala  Swerling  Case  I)  or  Rayleigh  plus  an  constant  component 
of  equal  power  (ala  Swerling  Case  III),  and  by  its  diversity  order,  Ng,  defined  as  the  number 
of  independent  target  values  within  the  N samples  noncoherently  integrated.  It  is  always 
necessary  to  run  SNRN,  the  calculation  for  a nonfluctuating  target,  after  entering  N,  PF, 
and  PD  and  before  running  SNRF,  the  calculation  for  any  fluctuating  target. 


Specific  user  instructions  follow: 


STEP 

INSTRUCTIONS 

INPUT 

DATA/UNITS 

— 

KEYS 

L. 

OUTPUT 

DATA/UNITS 

1 

Enter  program  SNRN 

2 

(Repeat  as  desired) 

N 

BSSalnH 

PD  > 0. 3 

ISTP]  VI  J 

PF  ^ 0.5 

(otqIUU 

La.1[7Z 

SNRN  dB 

If  desired 

r.-  'wsj 

Gi  dB 

To  include  collapsing  loss  or  tgt  fluct  go  to  step 

3 

l ::n  i _ j 

3 

Enter  program  SNRF 

r.  j izzi 

Samples  containing  signal 

1<NS-N 

LST'oJ  l_6_  J 

Order  of  target  diversity 

l<NesNs 

fsToj  n i 

Go  to  step  4,  5,  or  6 

LZ7  !ZZ 

4 

For  nonfluctuating  target 

laz  l:z 

SNRC  dB 

5 

For  Rayleigh  const  tgt  (Inc  Cases  III  & IV) 

[B  ] [_□ 

SNRF  dB 

6 

For  Rayleigh  & const  tgt  (Inc  Cases  III  & IV) 

rc  . 7 ! 

SNRF  dB 

If  desired 

CRTs  J [ZZ 

Lf  dB 

!R/s  1 1 j 

Gd  dB 

(Steps  4, 5,  and  6 may  be  repeated  in  any  order) 

R/Ro 

i 


j 

i 


I 


i 

1 

! 


i * 


Gd  : Diversity  gain  (dB) 

Gi  : Integration  gain  (dB) 

Lf  : Fluctuation  loss  (dB) 

N : Signal  and/or  noise  samples  integrated 
For  Cases  I and  III  : Ne=l 
For  Cases  n & IV  : Ne=Ns 
Ns  : Samples  within  N containing  signal 
PD  : Probability  of  detection 

PF  : Probability  of  false  alarm 

R/Ro  : Ratio  of  detection  range  to  that  for  which  SNRf  = 0 dB 

SNRC  ; SNR  per  sample  for  nonfluctuating  target  w/collapsing  loss  (dB) 

SNRF  : Avg  SNR  per  sample  for  fluctuating  target  (dB) 

SNRN  : SNR  per  sample  for  nonfluctuating  target  (dB) 


5-9 


HP-65  Program  Form 

HP-65  SNRN  Program  Listing 


Sv.  ITCH  TO  W PROM  PRESS  ij  lPRGM  ! TO  Cl  EAR  MEMORY 


KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

REGISTERS 

0 

00 

4 ^ 

04 

Rl  N 

ST08 

3308 

>T 

04 

RCL2 

3402 

PD 

tr  ..  j 

08 

b2  = 0.0448 

2 

02 

X 

71 

b2t 

R-)  PD 

X 

71 

• 

83 

1 

01 

9 

09 

gx>y 

3524 

1 > 2PD 

T 

U9 

3 

ST08 

3308 

2 

C2 

bi  = 0. 992 

£R1  . 

3508 

2PD 

+ 

. _ «1 

bl-d>2t  ...  

RCL2 
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PD 

X 

71 
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- 

51 
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1 

01 

xi 

B 

12 
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_+ 

61 

l+bit+b2t2=D 

SNR1 

g 

35 

.gRt  . 

3509 

t 

Rs  SNRN 

DSZ 

83 

• 

83 

CHS 

42 

2 .../ 

02 

_gNOP 

3501 

7 

07 

Re  . . _ 

>ST04 

3304 

-Q-1(PD) 

i 

. oi 

ai  = 0.  271 

RCL3 

3403 

PF 

X 

71 

. . --j 

B. 

12 

Q_I(PF) 

2 

02 

Rr  

RC IA 

3404 

83 

- 

51 

iQ”^(PF)+Q”1{PDl 
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03 
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41 
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x 

71 
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+ 

61 
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2 

..  02 
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3507 

D 

-r 

L 81 

l-f 

81 

N/D 
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ST04 

3304 

xi 
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Lr_- 

51 

Q-^P)  J 

RCL1 

3401 

N 

Lsto 

33 

4- 

81 

XN 

| 9 .. 

09 

! \8ELS 

' E 

15 

RTN 

24 

9-1vP)  _ . 

A ^ 

|ST05 

3305 

SNRN  ! 

, LBL 

23 

B qV) 

ST08 

3308 

E 

15 

X 

c 

RCL4 

3404 

XI 

_L 

41 

D j 

E 

F~16  j . . . 1 

! 

.41 

f SNRX  1 

| ST04 

3304 

SNR1 

_t  - 

41 

0 . 

I RCL5 

3405 

9 

09 

1 

R/S 

84 

"SNRN 

• 

83 

LBL 

23 

2 

02 

3 ... 

E 

12 

P 

+ 

61. 

X+0.92  , 

4 _ _ _ 

_I 

41 

X 

P -JL 

X(X+0. 92)  .. 

5 ..  _ 

X 

71 

P2 

f 

6 

31 

_\T 

09 

n/X(X+0.  92) 

7 . . 

JjN  . 

07 

In  P2  . 

+ 

61 

x W“  ( 

8 

CHS 

42 

2 

02 

9 ...  . 

f 

31 

4- 

81 

j 

j£. 

09 

t 

^ - 

31 
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41 

LOG 

08 
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41 

1 
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41 
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00 

2 

83 

X 

71 

Q 
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HP-65  Program  Form 

HP-65  SNRF  Program  Listing 


SW'TCM  TO  W/PflGM  PRESS  ' fj  [ PRGMJ  TO  CLEAR  MEMOR1 


KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

REGISTERS 

LBL 

23 

X 

71 

Ri  N 

A 

11 

2 

02 

0 

00 

- 

51 

D34 

GTO 

22 

LBL 

23 

R?  PD 

1 

01 

* 

2 

02 

D 1 

LBL  ! 

23 

f „4 

.31  J 

B 

12 

LOG 

08 

R3  PF 

RCL3 

3403 

PF 

1 

01 

f “'1 

_ 31 

0 

. . 00 

LN  1 

07 

X 

71 

10  log  D ...  

R4  SNR1 

RCL2 

3402 

PD 

RCL4 

3404 

SNR1 

f 

31 

- 

" 51  H 

LN 

07 

RCL7 

3407 

Ne 

Rl,  SNRN 

4- 

81 

In  PF/In  PD 

-r 

81 

Lf 

SNRF  1 

1 

01 

LBL 

23 

1 

- 

51 

D12 

- 

!i  . 

01 

Lf 

He  Ns  ] 

GTO 

22 

It 

41 

■) 

2 

02 

h 

,.j.  - - 

41 

" 

'LBL 

23. 

P.CL1 

3401 

N 

R7  Ne 

C 

13 

RCL6 

3406 

Ns 

— 

• 

83  _ 

•f 

81  j 

3 

03 

_ - _ . . . 

f 

31 

R8  SNRN 

6 

06 

'log 

08 

1 

01 

1 

01 

RCL3 

3403 

PF 

: o 

Q0. 

r> 

n9  _ 

f 

31 

1 

IX 

71 

io  log  (N/Ns) 

LOG 

08 

1 

1 

+ 

61 

- 

ui 

(.  361  - log  PFJ 

RCL8 

3408 

SNRN 

LABELS 

l5  _ ... 

03 

I 

+ 

61 

A Nonfluct 

• 

83 

ST05 

3305 

B ~Sw  1&2I 

2 

02 

— 

R/S 

84 

jSNRF 

C Sw  3&4 

7 

07 

3508 

D 

RCL2 

3402 

' PD 

R/s  . 

84 

Lf  J 

E 

CHS 

42 

JRCL7 

3407 

Ne 

O . . . 

1 

01 

1 

01 

1 Used 

+ 

61 

- 

. _ 51 

Ne-1 

? Used 

.31 

X 

71 

3 .... 

fj 

09 

i 

R/S 

84 

Gd 

4 

4- 

81 

3. 21/41- PD 

1 

-LBL 

23 

5 . 

CLSTX 

3500 

>4  1-pd 

1 

IE 

15 

6 

• 

83 

RCL5 

3405 

SNRF 

7 . 

9 

09 

CHS 

42 

8 

{6 

06 

4 

04 

9 

X 

71 

0 

00 

— 

51 

(3. 27A/1-PD)  - 0. 96*J1- 

•PD 

4- 

.fll 

FLAGS 

1 

01 

fi 

32 

1 

6 

83 

LOG 

08 

2 

02 

[R/S 

84 

R/Ro  _......  ... 

2 

9 

09 

1 

T) 

1 51 

u J 

PC 
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L 


( 


j 

i 


i 

i 

I 
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3,  CFAR  DETECTION,  RECURSIVE  SOLUTIONS 

These  HP-65  programs  calculate  the  detection  probability,  given  the  average  sample 
SNR  and  target  model  for  an  adaptive  detector  threshold  which  *s  set  proportional  to  the 
noncoherent  integration  of  R noise  samples.  These  programs  require  initial  calculation  of 
the  threshold  proportionality  constant,  T,  or  equivalently  A = 1/(T+1).  Two  cards  must  be 
entered  for  the  P2  case  or  three  cards  must  be  entered  for  the  general  case.  The  partition- 
ing is  such  that,  the  first  program  calculates  TO,  a starting  value  cf  an  iterative  solution  for 
T.  The  HP-65  P2C  program  does  the  iteration  for  T and  also  calculates  P2.  By  rerunning 
HP-65  TO  with  a given  P2,  the  required  average  SNR  can  also  be  calculated  by  HP-65  P2C. 

The  programs,  PGC(l)  and  PGC(2),  are  used  with  TO  to  calculate  detection  probability 
for  the  general  chi-squared  target  model. 


Specific  user  instructions  follow: 


STEP 


INSTRUCTIONS 


INPUT 

DATA/UNITS 


KEYS 


OUTPUT 

DATA/UNITS 


2 

3 


Enter  program  TO 


_N_ 

K 


PF 


Go  to  step  2 or  6 


00  z 

1 :n  z 
l a;:j  z 
zz 


□ 

□ 

\] 


Enter  program  P2C 


Go  to  step  3 or  4 


For  P2 


SNR  dB 


z 
ZZ1  z 
[ B..  j z 

z:j  z 

lzj  z 

iB/si  z 

Repeat  steps  4 & 5 for  new  P2  or  to  to  step  3 o|r  6 as  desifretj  1 V 


Repeat  for  new  SNR  or  go  to  step  4 or  6 as  desped 


For  SNR,  given  P2,  enter  program  TO 


P2 


Enter  program  P2C 


□ 

□ 


TO 


P2 


T20 


SNR  dB 


For  general  target  model  enter  program  PGC( 


SNR  dB 


K 


Enter  program  PGC(2) 


For  new  SNR  or  K enter  PGC(l)  and  go  to  step 


[R/sj  C 
CZJ  r 
dc  r 
mm  c 

zzz. 


A 

K 

N 

PF 

P 

PO 

P1-P4 

R 

SNR 

T 

TO 

T20 

X 


1/ (T+l) 

Chi-squared  distribution  parameter 

Signal  and/or  noise  samples  integrated 

Probability  of  false  alarm 

Probability  of  detection  for  chi-squared  target 

Probability  of  detection  for  nonfluctuating  target 

Probability  of  detection  for  Swerling  Cases  1-IV 

Noise  samples  integrated  to  set  threshold 

10  log  X 

Threshold  setting  di/.sor 
Iterative  solution  start  value  for  T 
Iterative  solution  start  value  for  T2=T/(1+X) 
Average  sample  S/N  within  N samples 


Zase  0 : K-105  ( 3 place  acc) 

Case  I : K=1 
CaseH  : K=N 

7 (Genl  Rayleigh  target  : K=Ne) 

\ Case  IH  : K=2 
Case  IV  : K=2N 

(Genl  Rayleigh  + equal  constant  target : K=2Ne) 
Weinstock  : 0 < K < 1 
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HP-65  Program  Form 

HP-65  TO  Program  Listing Page_ of 


S.-.irCH  TO  W/PfiGM  PRESS  I ] ' PRGM  TO  CLEAR  MEMORY 


KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS  j 

t 

KEY 

ENTR’i 

CODE 

SHOWN 

COMMENTS 

REGISTERS 

JLfiL  .. 

^ i . 

1 P2 

2 

Ri  N 

C 

l C 

P2 

f 

+ 

R+.922,  R-l 

SF1 

i 

T 

B 

R?  * 

.£T04L. 

P2 

X 

BLL.  ...  __ 

RCL3  ! 

TF  { SAVE  TF 

i 

ST08 

I IN  R 5 

. _j 

rgLSTx 

B,  1,B(  ) 

R,  Tor 

GTO 

_ J 

- 

hr 

; 1 

rclT4 

N 

LBL 

N 

X 

(1-B)  N,  B(  ) 

r4  hn 

A 

.. . 4 . 

t .. 

! + 

" 

11.  . 

ST04 

PF,  R,  N 1 

® 

1 1 

S 1 j 

t 

RCL4 

.. 

PIN 

r5  . . 

STG2 

R 

PF 

RCL2 

R 

41  ..1 

A 

S — 

STOl 

N 

i 

. i 

l 1/X  . 

Re 

LBL 

! 

8 

pinI/B,  1,4 4 'J 

ur.": 

i 

; yx  , 

RCL4 

PIN 

i 

. - 

R7 

it  . 

jgLSTx 

LOG  ! 

i 

4_ 

/ 1 -PIN1/*1)/ PIN1/!!,  [ J 

[CHS 

L 

X 

1/TO 

Rs  tF 

[ RCLl 

_ N. 

1 

.g  . - 

h t 

i 

.1/X 

; f 

ST03 

TO 

Rq 

n/n 

R/S 

s - -j 

L 

N-n/N.  L 

i 

gLSTx 

r 

NN 

- .j 

. 

LAOELS 

1 

j 

. 

A PIN  — PF 

’ *. 

n/n-i.  N-n/N.  I 

1 

B J 

IgRt 

L 

r.PiN— P2 

ml  _ 

1 

D 

^ r 

nTl 

i 

C 

. 

»/!  + */!_-  i 

-! 

0 

_gLSTx 

n/l 

i 

1 . . 

(vL  + vN-  1/vL,  NWN.  Li 

) 

2 . 

_gx— y 

j 

3 1 

gRi 

L.  N-n/n 

4 

-r 

5 . 

2 

| 

6 

. 

j 

7 . 

A 

J 

' 

fl 

(N-\fa)/2. 3L.  In/lWn-1)/'/ 

ti 

9 

i + 

j RCL2 

R 

FLAGS 

1 

| ; 

1P1N-P2 

- 

R-l 

RCL2 

R 

2 

• 

8 

rn, 
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HP-65  Program  Form 
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SWITCH  TO  W'PRGM  PRESS  f 1 PR(.M  ] TO  CLEAR  MEMORY 


KEY 

ENTRY 

CODE 

SHOWN 

C OMMENTS 

KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

REGISTERS 

RCL3 

TO 

f1 

Ri  N 

D 

P 

LOG 

RCL4 

PIN 

1 

•f 

P/PIN 

+ 

R?  R 

f 

RCL3 

T 

LN 

In  (P/PIN) 

X 

T2-T(l+X) 

RCL6 

0 

D 

Rs  T 

RCL5 

P 

R/S 

P2 

-r 

O/P 

LBL 

KCL2 

R 

JJ 

T _ 

— 

£ 

QC 

RCL3 

T 

J _ . . . . 

~ 

R/T,  Q/P,  ln(P/PIN) 

- 

1 .. 

_J 

Rs  PMS 

T 

AT 

+ 

STO 

IST07  ; 

T+l 

+ 

T* 

J J 

Re  QMS 

3 

(T-T+AT) 

RCL2  \ 

R 

RCL3 

T 

Ls  . . 

~ 

yx 

PM—  f T/(T+i)l R 

R7  T+l 

fL 

ST05 

PMS-PM 

ABS 

0 

M— 0 

EEX 

STOS 

QMS-0 

Rb  T 

CHS 

LBL.  . _ 

6 

10-».  IAT/TI 

J __  j.._  „ 

M,  PM  ...  

gx.js  y 

RCL2  ! 

R 

Rq  lOgiC 

GTO 

i&~~y  ! 

M,  R,  PM 

0 

RETURN  TO  PRGM  START 

+ 

rR+M 

RCL3 

T 

i 

gLSTx 

M 

LABELS 

r1 

i 

iel  J 

R+M,PM,  ~ M_  . ..  ... 

A . 

TFl 

!X 

. . . _J 

B SNR— P2 

R/S 

u. 

RCL7 

T+l 

C 1 

gNOP 

M-  PM— PM  (R+M)/ (T+ll 

D T—P 

r1 

STO 

E 

SF1 

+ 

nPRGMST 

RCL8 

T,T2  ^ RESTORE  T 

6 

(QMS— QMS  + M-PM) 

1 

ST03 

) IN  R3 

gRt 

M 

2 

4- 

1 

3 

i 

+ 

M — M+l 

4 

- 

X 

RCL1 

N,  M,  M-  PM 

5 

f 

gx— y. 

fi 

LOG 

RCL5 

7 

1 

RTN 

P-PMS 

a 

0 

. 

J&J&lPM  

0 

X 

PM 

a/s 

SNRdB 

STO 

FLAGS 

LBL 

+ 

1 P2—ITN 

B 

SNRdB 

5 

/PMS-PMS  + PM) 

1 

gLSTx 

M.  PM 

2 

0 

1 

GTO 

i-i  ziz.: 

1 
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HP-65  Program  Form 
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SWITCH  TO  W'PRCM  PRES 3 “ PRGM  ] TO  CLEAR  MEMORY 


KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

KEY 

ENTRY 

CODE 

SHOWN 

COMMENTS 

REGISTERS 

LBL 

RCL5 

P-PMS 

Ri  N 

0 

Q/P 

RCL3 

RCL2 

R 

RCL3 

RCL3 

T 

R?  R 

1 

i.t.t 

-r 

R7t,  Q/P,  In  (P/PIN) 

+ 

- 

ST07 

T+l 

T 

AT 

R^  T 

-r 

STO 

A 

RCL2 

R 

+ 

K 

3 . . 

/T—  T+AT} 

r4  pin 

PM—  [ T/  (T+l)  ]R 

RCL3 

T 

Z 

ST05 

PMS—  PM  1 

^ “1 

0 

M— 0 

! g J.. 

Rs  PMS 

ST06 

QMS-0 

ABS 

| ^T/T  | 

K 

LBL 

EEx 

1 

M,  PM 

• CHS 

Rb  QMS 

RCL2 

R 

ilO"6,  |AT/T| 

M 

gx—y 

r M,  R,  PM 

gxs  y 

1 + 

R+M 

GTO 

Fw  T+r 

[gLSTx  T 

M 1 

0 

| 

PM 

gRl 

R+M,  PM,  - , M 

RCL3 

T 

hr 

1 

RS  PMS 

RCL7 

1 T+l 

+ 

» 

• 

| M-PM 

g 

STO 

1/X 

R9  logic 

+ 



ST03 

i A-1/(T+1) 

6 

(QMS— QMS+M*  PM) 

R/S 

SNR 

IgRt 

M 

LBL 

t 

1 ABELS 

1 

A ~ 

— J J 

K 

ASNR,Kenti 

* 

M-M+l 

ST  05 

K ^ 

A 

B ... 

RCL1 

1 N,  M i 

gRl 

SNR 

C .. 

1 

D 

. 

GTO 

r 

0 

F 

2 

gRl 



M.M-PM 

— 

— 

— 

0 . 

1 

i PM 

_LOG  _ 

X 

? 

STO 

i __  . _____ 

RCL1 

J 

3 

+ 

X 

4 

5 

(PMS— PMS+PM) 

ST04 

z 

5 

gLSTx 

M,  PM 

1 

6 

GTO 

RCL3 

A ......  . . 

7 „ 

1 

- 

R 

LBL 

RCL2 

R 

9 

2 

L£L  . .. 

RCL5 

P—  PMS  1 

5* 

FLAGS 

RCL4 

' PIN  ! 

ST07 

' PM-d-A)^ 

1 

~ 

ST08 

PMS- PM 

f 

1 

7 

' LN 

r In  (P/PIN) 

ST06 

1. 

M— 1 

IRCL6 

Q— QMS 

,R/S 

1 
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SECTION  VI 
HP-67  PROGRAMS 

The  HP-67  is  very  similar  to  its  predecessor,  the  HP-65.  As  a result,  programs 
written  for  the  HP-65,  such  as  theprograms  of  this  report,  can  almost  be  transcribed  one- 
to-one  for  the  HP-67.  The  major  difference  is  the  greater  memory  and  programming 
capacity  of  the  67.  This  was  used  to  combine  and  store  multiple,  related,  HP-65  programs 
on  single  HP-67  program  cards. 

The  programming  differences  between  the  two  calculators  that  prevent  exact  one-to- 
one  transcription  are  noted  here  for  future  reference: 

• MERGED  INSTRUCTIONS 

The  HP-67  has  more  merged  instructions;  e.g. , "STO  + 8"  on  the  HP-67 
requires  3 program  lines  on  the  65. 

• CONDITIONAL  BRANCHING 

The  HP-65  skips  over  two  program  steps  if  the  conditional  test  is  false.  The 
HP-67  only  skips  one  program  step  when  the  test  is  false. 

• INDEX  REGISTER 

The  HP-65  uses  register  R6  as  an  index  register.  The  index  register  in  the 
HP-67  is  denoted  "I".  It  can  be  used  for  real  number  storage  as  well.  The 
register  R8  in  the  HP-67  is  for  data  storage  only. 

• PROGRAM  STORAGE 

The  HP-65  has  a capacity  of  100  program  steps.  The  HP-67  has  224. 

• DATA  STORAGE 

The  HP-65  has  9 storage  registers  including  register  9 which  is  net  fully 
available  because  it  is  used  for  internal  subroutines.  The  HP-67  has  26  data 
storage  registers  including  the  "I"  register.  All  are  fully  available. 

• LABELS 

The  HP-65  has  15  labels  for  program  entry  points,  subroutines,  and  branch 
points.  The  HP-67  has  20. 

The  HP-67  detection  programs  are  not  as  completely  annotated  as  the  HP-65  pro- 
grams. However,  the  HP-67  programs  can  be  readily  related  to  the  corresponding  HP-65 
programs.  Program  steps  and  labels  that  are  different  are  noted. 


1. 


FIXED -THRESHOLD,  RECURSIVE  PROGRAMS 


These  HP-67  programs  calculate  the  detection  probability,  given  the  number  of 
samples  noncoherently  integrated,  the  false-alarm  probability,  and  the  average  sample 
signal-to-noise  ratio  (SNR)  (in  dB)  for  the  various  target  models.  Unlike  the  HP-65  pro- 
grams, each  program  includes  the  threshold  determination.  The  PI  and  P2  calculations 
are  also  included  on  a single  program.  In  addition,  the  P1-P2  program  can  also  calculate 
by  iteration  the  value  of  SNR  required  for  a Case  II  target  with  a given  P2. 

Specific  user  instructions  follow: 


STEP 

INSTRUCTIONS 

INPUT 

DATA/UNITS 

KEYS 

OUTPUT 

DATA1UNITS 

EB 

Enter  PO,  P1-P2,  P3,  P4  or  PG  program 

N 

i “i  r:zi 

l a . ] | ;;zi 

[ i czj 

PF 

PF(calc) 

Go  to  step  2,  3,  4 or  5 as  appropriate 

2. 

For  PO,  P2,  P3  or  P4  do  either: 

SNRdB 

r b_..i  L_n 
i rj  ra 
i - "I  tzz 

P 

or: 

X 

P 

Repeat  as  desired 

For  P2  go  to  step  3 if  desired 

1"  icn 

3. 

For  SNR  or  X given  P2  (P2  prgm) 

P2 

■tag  — 

Hi  HI 

SNRdB 

X 

Repeat  or  go  to  step  2 as  desired 

EB 

For  PI  do  either: 

SNRdB 

■ran 

Hil 

MBIBU 

PI 

or: 

X 

PI 

Repeat  a3  desired 

5. 

For  P6  do  either: 

SNR  dB 

i b ;i  zz 
[T-J  CEl 

K 

or: 

X 

K 

6. 

Enter  desired  K if  different  from  display 

K 

rzj  r~j 

P 

Repeat  or  go  to  step  5 as  desired 

i 

i 

■ 

After  steps  2.  3.  4 or  6 do  anv  of 

rzi  lzj 

the  following  if  desired: 

(RCL]  L 0 ] 

PF  (calc) 

[RCL]  LX1 

[rcl]  r~r\ 

N 

(May  be  done  after  step  1) 

Y 

_ . 

|RCL|  | 9 | 

(RULJCT 

(Rcl]  rrn 

iRCL]  [ u 

P 

X 

SNRdB 

(PG  program  only) 

K 

K : Chi-squared  distribution  pai  aneter 

N : Signal  and/or  noise  samples  integrated 
Ne  : Target  diversity  within  N san  pies 
PF  : Probability  of  false  alarm 

P : Probability  of  detection  for  chi-squared  target 

PO  : Probability  of  detection  for  nonfluctuating  target 

P1-P4  : Probability  of  detection  for  Cases  I-IV 
SNR  : 10  log  X 

X : Average  sample  power  S/N  within  N 

Y : Fixed  detection  threshold 

Y2  : Y/(l+X) 


Case  0 : K=105  (-  3 place  acc) 

Case  I : K=1 

Case  n : K=N 

(Genl  Rayleigh  target  ; K=Ne) 

Case  HI  : K=2 
Case  IV  : K=2N 

(Genl  Rayleigh  + equal  constant  target. K=Ne) 
Weinstock  : 0 < K < 1 
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HP-67  PO  Program  Listing  (Cont) 


STEP  KEY  ENTRY  KEY  CODE 


COMMENTS 


STEP  KEY  ENTRY  KEY  CODE 


COMMENTS 


HP-67  P1-P2  Program  Listing 

STEP  KEY  ENTRY  KEY  CODE  COMMENTS 


LBLA 


CHS 


STOl 


LSTX 


+ 


LSTX 


X 


ST02 


LB  LI 


E 


RA 


RCL4 


X 


STO  + 2 


GTOl 


RCL8 


F ? 2 


GTOO 


STOO 


RTN 


LBLO 


ST09 


71 


61 


3302 


312501 


312215 


3408 


3554 


81 


3271 


2201 


3408 


357102 


2200 


3300 


3522 


312500 


312515 


02 


STEP  KEY  ENTRY  KEY  CODE 


RCL 


X s Y? 


RCL2 


Ra 


X 


RA 


STO  + 8 


R 


GTCK6 


LBLB 


0 


STOA 


LOG 


0 


STOB 


mm 


X— Y 


ST02 


RCL6 


ST02 


RCL8 


ST  09 


RTN 


LBLC 


SF2 


RCL2 


ST 


RCL1 


m 


2224 


8 


312512 


3253 


322512 


3311 


3153 


01 


00 


3309 


3522 


312513 


3 


312211 


COMMENTS 


REGISTERS 


1 N 

2 Y 

3 Z 

■a 

5 

65fTEMP 

| 

SI 

S2 

S3 

S4 

S5 

S6 

S7 

SNRdB 
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I 


HP-67  P1-P2  Program  Listing  (Cent) 


ST'iP  KEV  ENTRY  KEY  CODE 


COMMENTS 


STEP  KEY  ENTRY  KEY  CODE 


COMMENTS 


HP-67  P3  Program  Listing 

KEY  ENTRY  KEY  CODE  COMMENTS 


STEP  KEY  ENTRY  KEY  CODE 

rool  [ LBL 


STEP  KEY  ENTRY  KEY  CODE 


61 


COMMENTS 


-JP-67  P3  Program  Listing  (Cont) 


STEP  KEY  ENTRY  KEY  CODE 


COMMENTS 


STEP  KEY  ENTRY  KEY  CODE 


COMMENTS 


HP-67  P4  Program  Listing 


STEP  KEY  ENTRY  KEY  CODE 


COMMENTS 


LBLA 
-STD4 
LOG 
CHS 
X-—  Y 


312511 

3304 

3153 

42 

3552 


STEP  KEY  ENTRY 

zzz_r.cx.i~: 

X^Y  ? 

50  RTN 
RCL2 

ILC 

X 

r-  ^ A 


KEY  CODE 
61 


COMMENTS 


HP-67  P4  Program  Listing  (Cont) 

STEP  KEY  ENTRY  KEY  CODE  COMMENTS 


X 


STEP  KEY  ENTRY  KEY  CODE 


COMMENTS 


X 


STO  + 


CLX 


RCL1 


2 


LSTX 


RCL4 


RCL3 


ST04 


+ 


RCL6 


+ 


ST  06 


RCL1 


X 


GT04 


RCL5 

3405  ^ 

GT03 

2203 

ITTI  — 

'512504 

RCL8 

IHEOiniH 

m ■ 

3309 

RTN 

3522 

****tmm 

mymwm 

RCL6 

3406 

RCL1 


160  I X ^ Y ? 


j A~N  1 PF— • [fiSNRdB— P |c 


LABELS 

Id 


FLAGS 


oop  Branch 


SET  STATUS 

FLAGS 

TRIG 

DISP 

ON  OFF 

0 □ □ 

DEG  □ 1 

1 FIX  □ 

1 □ □ 

GRAD  □ 

SCI  □ 

2 □ □ 

RAD  □ 

ENG  □ 

3 □ □ 

n 

k 


HP-67  PG  Program  Listing 

STEP  KEY  ENTRY  KEY  CODE  COMMENTS 


001 


CHS 

42 

xsy 

3552 

STOl 

3301 

\wm 


LSTX 


X 


2 


3 


X 


+ 


2 


LBL1 


E 


STO  + 2 


LSTX 


ABS 


EEX 


RCL2 


3302 


312501 


312215 


408 


3554 


336102 


6 

06 

XSY? 

3271 

GTOl 

2201 

RCL8 

3408 

STOO 

3300 

12515 


01 


imm? 

m i 


PFC  N c Y 


KEY  ENTRY 

KEY  CODE 

\ + 

1 wn 

rsh 

— 

RTN 

3522 

RCL2 

3402 

RA 

3554 

X 

71 

RA 

3554 

STO  + 8 i 336108 


Ra 


GTOfi 


R 


LBLB  I 312512 


01 


10* 


LBLb 


81 


3253 


322512 


0 

00 

X 

71 

STOB 

3312 

RCLA 

3411 

RCL1 

3401 

STG13 


RCLE 


S 


RCL2 


C 


ST  07 


ST06 


LBLO 


RCL1 


LBL2 


RCL8 


STOS 


3313 


3415 


84 


312513 


3315 


3402 


3 


01 


3306 


12500 


REGISTERS 

HPF^ls 

SU1 


S4  S5 


■Mill 


+ 


STOS 


6 M 


IS6 


312502 


3 


3305 


01 


13 


3413 


WiO  ! 

in 


COMMENTS 


8 SNRdB 

Lc  ,z 

D 

E K 
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11^-67  PG  Program  Listing  (Cont) 

STEP  KEY  ENTRY  KEY  CODE  COMMENTS 


51 


IlilSUL 


ST04 


LSTX 


Rv 


RA 


LBL3 


GSBa 


RCL1 


RCL4 


RCL3 


X 


Hi! 

iwmm 


GT03 


■ ■ 

LBL4 

r i 

RCL5 

ST07 


51 


3554 


312503 


322211 


2203 


312504 


3405 


STEP  KEY  ENTRY  KEY  CODE 


1 I 01 


ST06 


RTN 


COMMENTS 


LABELS 


*N  PF-~  %NR  d£—  c K— ■ P 


d sub  r x-  r 


ftl"  LOOP!1  Y LOOP  BRANCH 


FLAGS 


Branch 


! SET  STATUS  i 

FLAGS 

TRIG 

DISP 

ON  OFF 
ODD 

DEG  D 

FIX  D 

i D D 

GRAD  D 

SCI  D 

2 D D 

RAD  □ 

ENG  D 

3 D D 

n 

2.  FIXED-THRESHOLD,  BARTON  ALGORITHM 

The  program,  HP-67  SNR,  calculates  the  required  average  SNR  in  dB  for  a given 
detection  probability  and  target  model.  The  target  model  is  specified  by  its  probability 
density  function,  i.e. , nonfluctuatimg,  Rayleigh  (i.  e. , Swerling  Case  I)  or  Rayleigh  plus 
an  equal  power  constant  component  (i.e. , Swerling  Case  II)  and  by  its  diversity  order,  Ne, 
defined  as  the  number  of  independent  target  values  within  the  N samples  noncoherently 
integrated. 

Specific  user  instructions  follow: 


STEP 

- - 1 --  - 

INSTRUCTIONS 

iNPUT 

DATA/UNITS 

KEYS 

OUTPUT 

DATA/UNITS 

1 

Enter  SNR  program 

2 

(Repeat  as  desired) 

N 

r_  ri  r:  i 

PD  s-  0. 3 

i:  r i ( i 

PFS0.5 

Lb'JLZ'.J 

SNRNdB 

For  integration  gain  if  desired: 

ircj  r.o  j 

Gi  dB 

3 

To  include  collapsing  loss  or  target  fluctuation 

I^NqSN 

rxi  rv.  \ 

]^Np-<N5 

r f j \d  ] 

Go  to  step  4.  5 or  6 

r 

4 

For  nonfluctuating  target 

LA..1I  _ 1 

SNRCdB 

5 

For  Ravleigh  target  tine.  Cases  I and  ID 

i..B.  11-  i 

SNRFdB 

6 

For  Rayleigh  constant  target  (inc.  Cases  HI  &IV 

i c j r i 

SNRFdB 

(Steps  4.  5 and  6 may  be  repeated  in  any  desire 

d order) 

i n.::i 

After  running  each,  if  desired; 

1R/S]  [ ..  . ] 

Ls  dB 

... 

IR/SJ  ( ...  J 

Gd  dB 

rn  r i 

R/b„ 

Go  : Diversity  gain  (dB) 

Gi  : Integration  gain  (dB) 

Lf  : Fluctuation  loss  (dB) 

N : Signal  and/or  noise  samples  integrated 
Ne  : Target  diversity  within  Ns  samples 
For  Cases  I and  HI  : Ne=l 
For  Cases  II  and  IV  : Ne=Ns 
Ns  : Samples  within  N containing  signal 
PD  : Probability  of  detection 

PF  : Probability  of  false  alarm 

R/Ro  : Ratio  of  detected  range  to  that  for  which  SNRF  = 0 dB 
SNRC  : SNR  sample  for  nonfluctuating  target  w/collapsing  loss  (dB) 
SNRF  : Average  SNR  per  sample  for  fluctuating  target  (dB) 

SNRN  : SNR  per  sample  for  nonfluctuating  target  (dB) 
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V 


HP-67  SNR  Program  Listing 


STSP  KEY  ENTRY  KEY  COPE  COMMENTS  STEP  KEY  ENTRY  KEY  CODE  COMMENTS 


6-14 


HP-67  SNR  Program  Listing  (Cont) 

^TEP  KEY  ENTRY  KEY  CODE  COMMENTS 


STEP  KEY  ENTRY  KEY  CODE 


COMMENTS 


‘A  SNRCdB- 


KB»8WBiataaMiawaiig^mg 


LABELS 

Hfc” 


FLAGS 


Subroutine  

1 Branch  2 Branch 


t Me  Subroutine 

4 2 

9 Branch  3 


SET  STATUS 

FLAGS 

TRIG 

DISP 

ODD 

DEG  □ 

FIX  □ 

1 □ □ 

GRAO  □ 

SCI  □ 

2 □ □ 

RAD  □ 

ENG  □ 

3 □ □ 

n 

i 

6-15/6-16 

3.  C FAR  DETECTION,  RECURSIVE  SOLUTIONS 

These  HP-67  programs  calculate  the  detection  probability,  given  the  average  sample 
SNR  and  target  model  for  an  adaptive  detector  threshold  which  is  set.  proportional  to  the  non- 
coherent integration  of  R noise  samples.  They  require  an  initial  calculation  of  the  threshold 
proportionality  constant,  T,  or  equivalently  A = 1/(T+1).  Each  program  includes  the  re- 
quired T or  A iterative  calculation.  The  HP-67  P2C  program  can  also  calculate  iteratively 
the  required  average  SNR  for  a Case  II  target  and  a given  detection  probability,.  The  HP-67 
PGC  program  calculates  detection  probability  for  the  general  chi-squared  target  model. 


Specific  user  instructions  follow: 


STEP 

INSTRUCTIONS 

INPUT 

DAT  A/UNITS 

Oi 

Enter  P2C  or  PGC  program 

N 

Fcr  P2C  go  to  step  2 or  3.  For  PGC  go  to  ste 


or: 

X 

Repeat  for  new  SNR  or  X or  go  to  step  3 as  des 

ired. 

3 

For  SNR  or  X given  P2 

P2 

Repeat  or  go  to  step  2 as  desired 


For  PG  SNR  dB 


Enter  desired  K if  different  from  display  K 


Repeat  4 as  desired. 


After  step  2,  3 or  4 do  any  of  the  following  desired: 


Note:  PGC  does  not,  store  P,  X,  or  SNR 
to  completion.  They  may  be  recalculated 
after  running  as  X=Z/N,  SNR  = 10  log  X 


1/(T+1) 

Chi-squared  distribution  parameter 

Signal  and/or  noise  samples  integrated 
Probability  of  false  alarm 
Probability  of  detection  for  chi-squared  target 
Probability  of  detection  for  nonfluctuating  target 
Probability  of  detection  foi  Swerling  Cases  1-IV 
Noise  samples  integrated  to  set  threshold 
10  log  X 

Threshold  setting  divisor 
Iterative  solution  start  value  for  T 
Iterative  solution  start  value  for  T2=T/ (1+X) 
Average  sample  S/N  within  N samples 


PGC  only  ( 


l A ... 

Ha  I 
LI  J 
(E£L\  Gl  ] 
ErcD  L4_  1 
[RCL]  Le  1 


p 


x 


SNR  (dB) 

A 

Z 

K 


Case  0 : K-10  (-  3 place  acc) 

Case  I : K=1 

Case  II  : K=N 

(Genl  Rayleigh  target  : K=Ne*) 

Case  III  : K=2 
Case  IV  :•  K=2N 

(Genl  Rayleigh  + equal  constant  target  : K=2Ne«) 
Weinstock  : 0 < K < 1 


*Ne:  Target  diversity  within  N samples 


HP-67  P2C  Program  Listing 


STEP  KE  eNTRY  KEY  COK 


001 


n sag  wmmmi 


COMMENTS  STEP  KEY  ENTRY  KT/  CODE 


X I 71 


COMMENTS 


RV 


ST02 


RV 


STOl 


LBL1 


RCL4 


LOG 


CHS 


RCL1 


LSTX 


3302 


3553 


3301 


312501 


3404 


3 


42 


3401 


41 


D 


RCL5 


LN 


STOC 


RCL6 


RCL 


01 


09 


42 


3533 


312500 


3403 


312214 


3405 


imm 


RCL2 


RCL3 


STO  + 3 


RCLC 


ABS 


EEX 


CHS 


6 


X s Y ? 


GTOO 


RCL5 


STOO 


TF3 


ST09 


RTN 


336103 


3413 


356 


43 


06 


3271 


22 


3405 


357102 


3300 


357103 


3309 


312512 


01 


0 


81 


LBLb 


RCL1 


X 


+ 


RCL4 


RCL2 


0 


X 


STOB 


RCLA 


LSTX 


0 

PF 

1 N 

2 

R 

3 

T 

so 

Si 

S2 

S3 

REGISTERS 


312214 


340 


6 SNRdB 

c 

D 

E 

|7 

T+l 

Sr 

xEMP 

9 

P2 

S7 

1 

S8 

S9 

L 


i 


| 


i 

i 

! 


) 


\ 

\ 

! 


HP-67  P2C  Program  Listing  (Cont) 


STEP 

KEV  ENTRY 

KEY  CODE 

COMMENTS 

STEP 

KEY  ENTRY 

KEY  CODE 

COMMENTS 

- 

STOS 

3309 

RTN 

3522 

170 

“ 

LBLU 

312513 

“ 

SF3 

355103 

ST04 

3304 

" 

RCL3 

3403 

“ 

ST08 

3308 

" 

*20 

GSB1 

312201 

“ 

JRCL3 

3403 

” 

RCL8 

3408 

STOS 

3303 

T- 

81 

190 

1 

01 

- 

51 

sioa 

3311 

LOG 

3153 



1 

01 

130 

0 

00 

X 

71 

_ ..  

STOB 

3312 

- 

RTN 

3522 

LBLD 

312514 

I'.'O 

I 

41 

i 

41 

1 

01 

+ 

61__ 

ST  07 

3307 

'40 

RCL2 

81 

3402 

. 

_YX 

3563 

— 

ST05 

3305 

— 

0 

00 

200 

ST06 

3306 

RCL2 

3402... 

X — Y 

3552 

+ 

61 

150 

LSTX 

Rv 

3582 

3553 

-X 

RCL7 

71 

3407 

81 

STQ  + 6 

-336  1 06 

210 

RA 

3554 

1 

01 

4* 

61 

■RCLl 

3401 

Y? 

3271 

160 

RTN 

. 35.22 

Rv 

3553 

STO  + 5 

_fii 

336105 

LSTX 

GTO/ii 

3582 

2l25 

220 



' 

_ 

- 

I A REM  C 

3 

S\  R,  PF- 

lNRdb-P 

CP— SNRdB 

° SUB 

E 

\o 

r L8U& 

SET  STATUS 

FLAGS  TRIG 

DISP 

■ 

b X — 

P 

C 

d 

e 

i 

ON  OFF 
ODD 

DEG  □ 
GRAD  □ 

FIX  S 

SCI  □ 

_ c 

LOOP 

1 LOOP 

2 

3 

4 

2 

A 

i □ □ 

6 

7 

8 

9 

C 

■ 2 □ □ 

3 □ □ 

RAD  □ 

ENG.  □ 
n v 

■ew 
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HP-67  PGC  Program  Listing 


COMMENTS 


HP-67  PGC  Program  Listing  (Cont) 

KEY  ENTRY  KEY  CODE  COMMENTS 


2200 


3/X 


ST03 


0 


RTN 


STEP  KEY  ENTRY  KEY  CODE 


3 I 


3522 


01 


HP  %%4MHH 

[_  3553  1 

—'ll  i m 

a 

322211  \ 

Rv 

mmm 

RCL9 

3409 

RCL7 

3407 

COMMENTS 


X 


ST  04 


RCLE 


im 


+ 

61 

LSTX  1 

3582 

X 


RCL5 


X 


5 


STO-9 


ST08 


2 


GTOl 


RCL4 


3308 


01 


3306 


312501 


3401 


X >Y? 


Rv 


RTN 


LB  La 


322511 


RCL2  ! 


RCLE 


x 


ST  05 


RCL8 


Rv 


Rv 


RCL3 


RCL7 


X 


ST07 


( i LBL6  i 31ZDU 
|A  Nt  R— - !%NRdB+K4C 


312506 


b 


s 


p 


RANCH 


LABELS 

ID  „ 


ST06 


FLAGS 


! SET  STATUS  1 

FLAGS 

TRIG 

D1SP 

ON  OFF 
ODD 

DEG  D 

FIX  D 

1 □ D 

GRAD  D 

SCI  D 

2 D D 

RAD  D 

ENG  D 

3 D D 

n 
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SECTION  VII 


ILLUSTRATIVE  EXAMPLES 

1.  HP-65  FIXED-THRESHOLD,  RECURSIVE 


Program 

Step 

Input 

Key 

Output 

Time  (m:s) 

Y-P2 

1 

N=10 

PF=10~6 

t 

A 

Y=32»  710341 

0:28 

2 

x=\/r o 

B 

P2=0. 733987 

0:08 

3 

P2=0,  9 

A 

Y2-45.  221305 

0:45 

R/S 

X=4. 257794 

0:01 

PO 

4 

X=n/10 

B 

PO^O. 853317 

2:23 

PI 

4 

x=VTo 

B 

P1=0. 485543 

0:18 

P3 

4 

x=\ITo 

B 

P3=0. 569375 

0:20 

P4 

4 

X=n/Io 

B 

P4=0. 781789 

0:33 

P6 

7 

K=105 

t 

X-n/10 

B 

P0«0. 853297 

2:19 

7 

K=1 

t 

x=^!To 

B 

P1=0. 485543 

2:35 

7 

K=10 

1 

x=4To 

B 

P2=0. 733987 

2:31 

7 

K=2 

! 

x=\fio 

B 

P3=0, 569375 

2:34 

7 

K=20 

t 

X=<n/Io 

B 

P4=0.  781789 

2:27 
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2.  HP-65  FIXED-THRESHOLD,  BARTON 


Program 

§t£E. 

Input 

Key 

Output 

Time  (m:s) 

SNRN 

2 

N=10 

STOl 

PD=0. 75 

STO  2 

PF=10~^ 

ST0  3 

A 

SNRN -4. 3^ 

0:09 

- 

Gj=7. 89 

3 

N -10 

S 

STO  6 

N =10 
e 

STO  7 

4 

A 

SNRC=4. 34 

0:01 

5 

B 

SNRF=4. 79 

0:03 

6 

C 

SNRF-4.57 

0:04 

R/S 

Lf=0. 22 

0:01 

R/S 

Gd-2.  05 

0:01 

E 

R/R  =0. 769 

0:01 

HP-65  CFAR 

DETECTION 

Program 

Step 

Input 

Key 

Output 

Time  {m:8) 

TO 

i 

N=10 

t 

R=16 

1 

A 

PF^IO-6 

T0=0, 224789 

0:04 

P2C 

2 

R/S 

T~0.  234428 

0:41 

3 

SNR=10 

B 

P2=0.  8676  5 2 

0:13 

TO 

4 

P2=0. 9 

C 

T20=2.  207242 

0:04 

P3C 

5 

R/S 

SNR=10. 339693 

1:10 

PGC(l) 

6 

R/S 

A=0.  810092 

0:41* 

7 

SNR=10 

t 

K=10 

A 

i 

0:01 

PGC{2) 

8 

R/S 

P2=G.  8676  52 

9:57 

PGC(l) 

7 

SNR=10 

t 

K-l 

A 

1 

0:01 

PGC{2) 

8 

R/S 

Pl-C. 564214 

10:42 

*0:14  if  step  2 has  been  run 
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4. 


KP-67  FIXED-THRESHOLD,  RECURSIVE 


Program 

Step 

Input 

Key 

Output 

Time(m:s) 

Any  below 

1 

N=10 

t 

PF=10~6 

A 

PF=0.  000001 

0:32 

P1-P2 

2 

SNR=5 

B 

P2=0. 733987 

0:11 

3 

P2=0,  9 

C 

SNR=6.  29x847 

0:56 

4 

SNR=5 

D 

P1=0. 485543 

0:24 

PO 

2 

SNR=5 

B 

PC:0. 853317 

2:46 

P3 

2 

SNR=5 

B 

P3=0. 569375 

0:08 

P4 

2 

SNR=5 

B 

P4~0. 781789 

0:39 

PG 

5 

SNR =5 

B 

6 

K-105 

C 

POisO.  853298 

2:40 

6 

II 

C 

P1=0. 485543 

2:59 

6 

K=10 

C 

P2=0.  733987 

2:55 

6 

K=2 

C 

P3=0.  569375 

2:58 

6 

K=20 

C 

P4=0.  781789 

2:51 

HP-67  FIXED-THRESHOLD,  EaRTON 

Program 

Step 

Input 

Output 

Time(m:S) 

SNR 

2 

N=10 

t 

PD=0,  75 

t 

PF=10 

D 

SNRN=4. 34 

0:13 

3 

Ns=1° 

* 

i 

N =10 

0 

fd 

4 

A 

SNRC=4. 34 

0:02 

5 

B 

SNRF=4. 79 

0:05 

6 

C 

SNRF=4. 57 

0:06 

R/S 

Lf=0.  22 

0:01 

R/S 

Gd=2. 05 

0:01 

E 

R/R  =0. 769 

0:01 
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HP-67  CFA.R  DETECTION 


Program 

Step 

Input 

Ke^ 

Output 

Time(m:s) 

P2C  or  PGC 

1 

N=10 

R=16 

PF=10“6 

A 

PF=0.  000001 

1:01 

P2C 

2 

SNR=10 

3 

P2=0.  867652 

0:14 

3 

P2=0. 9 

C 

SNR=10. 339693 

1:15 

PGC 

4 

SNR=10 

B 

K-10 

R/S 

P2=0. 86 76 52 

10:24 

4 

SNR=10 

B 

K=1 

R/S 

P1=0. 564214 

11:11 
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Please  make  the  following  changes  to  T1S  R79EMH5: 


p.  2-5 
p.  3-1 

P.  3-1 

p.  5-1 

p.  5-1 

p.  5-2 

p.  5-3 

p.  5-4 

p.  5 "4 

p.  5-7 

p.  5-7/ 
5-8 

p.  5-9 

p.  5-10 

p.  5-13 
p.  6-2 
p.  6-2 
p.  6-5 
p.  7-1 
p.  7-1 
p.  7-2 
p.  7-4 


Change  reference  number  for  Shniclman  from  2 to  4. 


Equation  (3-1)  center  should  read 


Q“\PF)  ~ Q 


Change  reference  number  for  Cann  from  8 to  7. 


Step  5 output:  Change  "D"  to  "P". 

in  instruction  table  - step  7,  first  key  block  after  K should  contain  "t 
Step  25,  Comment  should  read:  2.3  n/l  (sfh  + *sTn  - 1),  N - \!  N 
Step  3,  change  "3302"  to  "3303". 


Step  4,  change  "3403"  to  "3402". 


Step  19,  change  "3302"  to  "3402". 

On  program  line  52,  change  comment  fiom  "XB  • V/B"  to  "XB  • V/(M-N)". 
Step  5,  change  "3303"  *o  "3304". 


lnstr.  step  5 should  read:  "For  Rayleigh  tgt.  (Inc  Cases  1 & II) 

Line  21  comment  should  read  [q’\pF)  - Q_1(PD)J  ; Comment  column 
line  17,  remove  minus  sign.  L 

5 

In  nomenclature  table:  Case  0:  h 10 
Step  5,  "P6"  should  be  "PG". 

Add  to  step  1,  1st  block  , insert  t . 

On  program  line  089,  change  "X"  to  "+". 

Change  last  item  under  Program  from  P6  to  PG. 

Line  1,  step  1,  change  " 1 " to  "STOl". 

Par.  2,  step  3,  add  SNRF  in  Program  Column. 

Step  1 in  Key  column,  add  t in  Key  Column  after  N = 10  and  R - 16. 


